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(DETACHABLE  SUMMARY) 


THE  ANALYSIS  OF  THE  EFFECTS  OF  FRAME  RESPONSE 
ON  BASEMENT  SHELTERS  IN  TALL  BUILDINGS 


This  report  presents  the  results  of  a  program  to  develop  a  theoretical  analysis  of 
the  effects  of  frame  response  on  basement  shelters  in  tall  buildings.  The  objective 
was  to  determine  the  effect  on  a  upgraded  basement  key  worker  shelter  of  the 
aboveground  portion  of  the  structure  being  subjected  to  a  blast  wave  (30  to  50  psi) 
that  would  destroy  the  building. 

This  program  investigated  both  steel  and  reinforced  concrete  frame  structures 
with  the  most  emphasis  on  poured-in-place  reinforced  concrete  beam,  slab,  and 
girder  type  framing  and  the  poured-in-place  flat-slab  and  flat-plate  type  of 
construction.  These  types  are  very  common  in  the  National  Shelter  Survey  inventory 
of  upgradable  structures. 

A  prediction  technique  was  developed  using  both  hand  and  computer  analysis. 
This  technique  was  tested  using  a  previously  explosively  demolished  15-story  cast-in¬ 
place  reinforced  structure,  the  Continental  Life  Building  in  Atlanta,  Georgia.  The 
results  of  this  analysis  indicated  that  the  upgraded  basement  would  have  survived 
even  though  the  aboveground  portion  of  the  structure  was  exposed  to  50  psi. 

Large  portions  of  the  debris  landed  on  the  shelter  roof  so  that  serious  questions 
remain  as  to  the  advisability  of  using  such  structures  as  shelters  because  of  the 
possible  problems  of  entrapment  of  the  shelterees.  Other  problems  that  need  to  be 
addressed  are  the  large  amounts  of  dust  created  in  the  collapse  process  and  the 
possibility  of  fire  in  the  debris  pile. 


It  also  should  be  noted  that  the  primary  work  in  this  program  was  on  reinforced 
concrete  structures,  where  the  punching  effects  of  falling  columns  do  not  appear  to 
be  a  problem.  Based  on  the  experiences  of  an  explosive  demolition  contractor, 
Controlled  Demolition,  Inc.,  this  appears  to  be  a  serious  problem  in  steel  framed 
structures.  It  is  recommended  that  future  work  in  this  area  thoroughly  investigate 
the  collapse  of  steel  frames  both  theoretically  and  in  conjunction  with  building 
demolitions.  Also,  as  was  noted  in  Section  4,  very  little  information  is  available  on 
the  load  imparted  to  a  frame  by  failing  walls  at  pressures  above  15  psi.  Tests  at  the 
higher  pressures  of  interest  (30  to  50  psi)  need  to  be  conducted. 
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Section  1 

INTRODUCTION 


Civil  Defense  planning  in  the  United  States  is  currently  based  on  a  policy 
termed  "Crisis  Relocation".  This  policy  presumes  that  a  period  of  crisis  buildup  in 
the  world  —  similar  to  the  1961  Cuban  and  the  more  recent  Middle  East  crises  — 
will  precede  any  future  war.  This  period  of  crisis  would  allow  time  (a  few  days  or 
weeks)  to  conduct  a  number  of  activities  to  protect  the  civilian  population  and 
industry  from  possible  attack.  These  activities  include:  evacuation  of  the  major 
portion  of  the  population  to  low-risk  areas  where  only  fallout  and  possibly  low  level 
blast  protection  would  be  required;  the  hardening  and  protection  of  critical  industries 
and  facilities;  and  the  provision  of  shelters  for  a  small  contingent  of  key  workers 
(who  will  remain  behind  to  maintain  vital  services  -  communications,  fire  protection, 
military  production,  etc.). 

This  project  was  concerned  with  the  last  of  these  activities,  the  provision  of 
shelters  for  key  workers.  As  stated  in  the  RFP,  the  objective  of  this  program  was  to 
"Develop  a  theoretical  analysis  of  the  damage  to  key  worker  basement  shelters  that 
may  result  when  the  buildings  in  which  they  are  located  are  subjected  to  a  nuclear 
weapons  blast  (peak  overpressure  30  -  50  psi)  so  that  the  buildings  are  destroyed  and 
the  integrity  of  the  shelters  may  be  threatened.  The  response  of  the  building  frame 
as  it  relates  to  the  integrity  of  the  shelter  and  the  survival  of  those  in  the  shelter 
should  be  explored  in  detail." 

Various  key  worker  shelter  concepts  have  been  explored,  including  the 
upgrading  of  basements  in  existing  structures  (for  example,  Refs.  1  through  6.) 
Work  has  also  been  done  on:  predicting  the  collapse  levels  of  floors  over  basements 
(for  example.  Refs.  2  and  7  through  10);  failure  predictions  and  full-scale  tests  of 
walls  (for  example,  Refs.  11  and  12);  the  flow  of  blast  waves  through  openings  (for 
example,  Refs.  13  through  21);  and  the  effects  of  fire  on  basement  shelters  (for 
example,  Refs.  22  through  24). 
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Until  now,  however,  very  little  work  has  been  done  on  the  combined  problem; 
i.e.,  what  is  the  effect  of  the  building  collapse  on  the  survivability  of  an  upgraded 
shelter  located  in  the  basement?  There  are  several  concerns.  For  example,  the 
blast  wave  might  not  shear  off  a  building  just  above  the  basement,  but  might  cause 
localized  basement  wall  failures  or  upgraded  basement  roof  slab  failures,  because  of 
either  the  frame  response  or  the  impact  of  debris.  Moreover,  at'  pressure  levels  of 
30  to  50  psi  most  structures  will  collapse,  and  much  of  the  structure  may  land  on  the 
floor  over  the  basement.  The  buildings  considered  in  this  study  (4  to  6  story  and  10 
to  12  story)  are  expected  to  provide  the  greatest  opportunity  for  massive 
superstructure  collapse  onto  the  floors  over  basement  space.  In  addition,  structures 
of  these  story  heights  often  appear  in  groups;  thus,  if  debris  from  a  structure  falls 
clear  of  its  basement,  surrounding  structures  would  be  likely  to  supply  an  equivalent 
mass  of  debris. 

Relevant  past  research  is  best  represented  by  the  SRI  reports  by  Wiehle,  Refs. 
9  and  10,  and  its  limitations  are  best  illustrated  by  the  following  statements 
extracted  from  these  reports.  From  Ref.  9, 

"The  strength  of  the  exterior  walls  is  important  in  calculating  the 
collapse  of  the  frame,  since,  for  a  given  overpressure  level,  the  blast 
loading  on  the  total  wall  area  can  be  much  more  severe  than  the  blast 
loading  on  the  frame  alone  plus  an  impulse  loading  from  a  frangible- 
type  wall. 

To  investigate  the  relative  strength  of  the  exterior  walls  and  frame  of 
a  building  would  require  a  comprehensive  computer  program  that 
includes  inelastic  response  under  dynamic  loading  as  well  as  realistic 
frame  collapse  mechanisms.  Since  such  a  program  was  not  available,  a 
computer  program  for  analyzing  the  elastic  and  inelastic  dynamic 
response  of  two-dimensional  structural  frames  was  used  .... 
Although  the  program  does  not  include  frame  collapse  mechanisms,  it 
was  felt  that  the  results  would  provide  a  basis  for  estimating  the 
possible  collapse  strength  of  a  building  frame  relative  to  the  strength 
of  the  exterior  walls." 
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and  from  Ref.  10, 


"The  collapse  of  the  floor  slab  over  basement  areas  is  an  important 
consideration  in  determining  the  survivors  in  nuclear  blast 
environments.  However,  collapse  predictions  for  the  floors  in  the 
Greensboro- High  Point  buildings  could  not  be  included  in  this  effort 
because  the  procedures  are  currently  being  developed.  The  analysis  of 
floor  slabs  will  be  included  in  the  building  collapse  predictions  when 
the  procedures  become  available." 

In  summary,  previous  efforts  have  consisted  of  investigations  to  establish 
exterior  wall  failure  overpressures  and  approximate  analyses  of  frame  response  in  a 
corresponding  range  of  overpressures,  but  excepting  the  SSI  work  on  the  key  worker 
manual  and  the  companion  technical  manual  (Refs.  1  and  2),  the  failure  analysis  and 
theoretical  efforts  have  been  performed  on  "as  is"  structures.  Such  analyses  fall 
short  of  present  needs  because,  as  discussed  in  Section  4,  the  failure  modes  are 
significantly  different  when  a  structure  is  upgraded.  Another  unfortunate 
circumstance  with  earlier  work  is  that  little  of  the  existing  test  data  on  wall/frame 
interaction  has  been  made  directly  available.  In  the  early  days  of  the  shock  tunnel 
test  program  (Ref.  11),  measurements  were  made  via  load  cells  installed  on  all  walls 
tested,  to  enable  this  interaction  to  be  investigated.  However,  objectives  at  that 
time  precluded  the  expenditure  to  reduce  the  unreported  data  for  current 
application,  so  that  only  reported  early  data  have  been  utilized. 

Recognizing  the  limitations  in  the  past  research  and  the  need  to  apply  the 
wall/frame  interaction  data  in  the  light  of  knowledge  of  basement  upgrading  schemes, 
the  research  described  in  this  report  was  performed  to  answer  the  following 
questions:  Given  a  time-varying  blast  load  on  a  frame  building  structure,  see  Figure 
1-1, 

o  What  are  the  various  conceivable  collapse  mechanisms  (e.g.,  "falling 
tree",  "soft  story"  and  fractured  members)? 

o  What  combinations  and  sequences  of  these  occur  and  which  are  the 
most  prevalent? 
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o  Given  the  prevailing  collapse  mode  sequence,  what  effects  do  the 
collapse  mechanism  and  debris  impact  have  on  the  upgraded  shelter 
space? 

The  technical  approach  taken  to  answer  these  questions  consisted  of  four  tasks: 
Task  1  -  Building  Selection 

Representative  medium-rise  and  high-rise  frame  buildings  with  basement  shelter 
space  were  reviewed  for  an  example  to  demonstrate  the  analysis. 

Task  2  -  Development  of  Upgrading  Plans 

Previous  SSI  shelter  work  was  used  to  prepare  shelter  area  upgrading  plans  to 
apply  in  the  structure  selected  for  the  example.  The  upgrading  was  for  combined 
blast  and  radiation  effects. 

Task  3  -  Theoretical  Analysis 

This  was  divided  into  three  sub  tasks: 

3A  -  Development  of  blast-generated  loading  criteria. 

3B  -  Development  of  joint  resistance  functions.  Ultimate  strength  capacities 
for  both  steel  and  reinforced  concrete  were  adopted  for  the  frame  mechanism 
joint  resistance  functions. 

3C  -  Theoretical  analysis  of  the  damaging  effects  of  frame  response.  The 
theory  of  plastic  analysis  was  used  to  predict  the  frame  failure  modes,  and  the 
resulting  force  and  distortion  information  was  used  to  evaluate  effects  on  the 
shelter  cover  slab. 

Task  4  -  Computer  Model  Development  and  Analysis 

The  dynamic  analysis  portion  of  the  program  used  the  static  and  dynamic 
features  of  the  Georgia  Tech  Integrated  Civil  Engineering  System  Structural  Design 
Language,  GTSTRUDL  (Ref.  25),  combined  with  a  progressive  collapse  mechanism 
model  to  predict  the  time  history  of  the  frame  response  and  mode  of  failure.  The 
results  were  then  evaluated  with  respect  to  the  effect  on  the  shelter  cover  slab. 
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The  remainder  of  this  report  is  organized  as  follows:  Section  2  discusses 
building  selection  (Task  1);  Section  3  deals  with  the  development  of  upgrading  plans 
(Task  2);  Section  4  presents  the  theoretical  analysis  (Task  3  and  its  associated 
subtasks);  Section  5  discusses  the  computer  model  development,  analysis  and  results 
(Task  4);  and  Section  6  presents  the  summary  and  conclusions.  Applicable  concepts 
of  plastic  analysis  are  discussed  in  Appendix  A,  computer  analysis  fundamentals  are 
presented  in  Appendix  B,  and  a  discussion  of  loading  functions  is  to  be  found  in 
Appendix  C. 


Section  2 

TASK  1  -  BUILDING  SELECTION 


INTRODUCTION 

The  selection  of  a  building  structure  for  analysis  required  that  each  candidate 
structure  be  evaluated  using  criteria  that  would  assure  that  the  final  selection  was 
representative,  compatible  with  all  of  the  tasks  required  in  the  program,  and 
ultimately,  would  lend  itself  to  a  definitive  assessment  of  the  computer  model.  The 
following  parameters  were  judged  to  be  the  most  important,  and  were  included  in  the 
selection  process: 

I.  The  structure  should  be  either  a  high-rise  (10  to  12  stories)  or  a 
medium-rise  (4  to  6  stories). 

II.  The  basement  should  be  large  enough  to  provide  shelter  space  for  50 
people  when  it  was  upgraded  in  accordance  with  Ref.  1  to  survive 
30  to  50  psi  —  before  giving  consideration  to  collapsing  stories 
above  the  floor  over  the  basement. 

III.  Sufficient  data  on  the  building  should  be  available  in  the  form  of 
structural  as-built  drawings  in  order  that  the  theoretical  analysis 
could  be  performed. 

IV.  The  structural  building  frame  should  not  just  "blow  away",  but  be  of 
a  type  that  offers  reasonable  lateral  resistance  so  as  to  develop 
representative  results  (e.g.,  a  prybar  effect  on  the  below-grade 
structure).  The  frame  columns  should  extend  through  the  ground 
floor  to  the  basement  foundation,  in  order  to  determine  the  failure 
effects  of  the  columns  on  the  ground  floor  slab. 
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V.  The  building  should  be  representative  of  those  in  use  today,  and  not 
be  a  unique  structure,  such  as  possibly  a  very  old  building  that  was 
designed  and  constructed  using  methods  and  materials  no  longer  in 
common  use. 

Two  concrete  structural  types  that  fit  these  criteria  are  shown  in  Figure  2-1. 


CANDIDATE  BUILDINGS 

Candidate  buildings  were  obtained  from  three  sources  and  a  preliminary  list  of 
structures  compiled  for  comparison  with  the  above  criteria  (items  I  -  V)  in  order  to 
select  one  for  complete  analysis.  The  three  pools  of  buildings  were:  (1) 
structures  from  the  National  Shelter  Survey,  (2)  structures  scheduled  to  be 
demolished,  and  (3)  structures  recently  demolished.  How  this  preliminary  list  was 
evaluated  is  discussed  below. 

National  Shelter  Survey 

It  was  believed  that  the  National  Shelter  Survey  would  be  a  viable  source  for 
candidates  for  analysis  because  of  the  large  number  of  buildings  and  the  various 
types  of  construction  that  were  included,  and  the  documentation  of  each  that  had 
previously  been  accomplished.  The  following  six  buildings  (see  Ref.  26)  from  the 
National  Shelter  Survey  were  selected  for  the  preliminary  list;  the  prints  and  data 
were  provided  by  James  E.  Beck  &  Associates. 

1.  Henry  R.  Landis  State  Hospital,  Philadelphia,  PA  (NSS  No.  110)  — 

Constructed  in  1960,  the  complex  consists  of  five  attached  structures,  three  seven 
stories  in  height,  one  three  stories  high,  and  one  two  stories  high.  The  basic 
construction  is  a  reinforced  concrete  frame  with  floors  and/or  roofs  consisting  of 
several  types  of  reinforced  concrete  construction,  including  one-  and  two-way  slab 
and  beam,  one-way  pan  and  joist,  waffle  slab,  and  open-web  steel  joists  with  metal 
deck  and  concrete  topping.  The  walls  are  of  masonry  construction,  generally 
concrete  block  with  brick  veneer.  One  of  the  seven-story  buildings  contains  a 
partial  basement. 
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These  structures  were  judged  to  be  unsuitable  for  this  program  because  of  the 
lack  of  size  (item  II). 

2.  Library  Building,  University  of  Illinois,  Chicago,  1L  (NSS  No.  149)  —  This 
building  is  a  four-story  structure  with  a  full  basement,  and  is  constructed  of 
reinforced  concrete  columns  and  girders  supporting  precast  concrete  beams  and  5-in. 
thick  slabs  concrete  topped  with  3-in.  of  concrete.  The  walls  are  of  infill  masonry 
construction.  The  building  was  constructed  in  1967. 

The  design  of  the  building  is  rather  unique  and  does  not  meet  the  criteria  of 
what  might  be  considered  "typical"  construction  (item  V).  The  building  frame  offers 
very  little  resistance  to  lateral  loads  (item  IV),  and  it  was  therefore  judged  that  the 
analytical  results  would  not  prove  satisfactory. 

3.  City  Hall,  Houston,  TX  (NSS  No.  196)  —  An  eleven-story  steel  frame 
structure  with  a  full  basement  constructed  in  three  stages  in  1937,  1939,  and  1953. 
The  photographs  of  the  building  and  the  NSS  Data  Collection  Form  indicated  that  the 
structure  was  a  possible  candidate  for  analysis;  however,  the  drawings  supplied  were 
completely  unreadable,  and  no  usable  data  for  analysis  could  be  extracted  (item  III). 

4.  Fidelity  Federal  Plaza  Office  Building,  Long  Beach,  CA  (NSS  No.  220) 

— This  structure  is  an  eleven-story  reinforced  concrete  frame  building  with  two  full 
levels  below  grade.  It  was  constructed  in  1967.  The  structural  floor  system 
consists  of  reinforced  concrete  one-way  slabs  and  beams,  and  the  walls  are  6-in. 
thick  reinforced  concrete.  The  building  is  very  large  and  complex  for  the  first 
efforts  at  collapse  prediction.  (This  is  not  one  of  the  criteria  for  rejection,  but  a 
somewhat  simpler  structure  was  sought.) 

5.  Eastland  Shopping  Center,  West  Covina,  CA  (NSS  No.  227)  —  This  shopping 
complex  consists  of  five  separate  structures;  i.e.,  one  three-story  department  store 
and  four  one-story  specialty  shop  buildings,  all  constructed  in  1956. 

The  department  store  is  of  reinforced  concrete  frame  construction  and  has  a 
small  partial  basement.  The  first  floor  is  of  reinforced  concrete  flat  slab  and  one 
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way  pan  and  joist  construction,  and  the  second  and  third  floors  we  all  one-way  pan 
and  joist  construction.  The  walls  are  a  mixture  of  reinforced  concrete  and  masonry. 

The  four  specialty  shop  buildings  are  all  of  the  same  construction  and  are 
located  over  one  full  basement.  The  first  floors  are  of  reinforced  concrete  one-way 
pan  and  joist,  one-way  slab  and  beam,  and  precast  prestressed  concrete  double-tee 
construction.  The  roof  is  entirely  constructed  of  double-tees,  and  the  walls  are  a 
combination  of  masonry  and  reinforced  concrete. 

The  department  store  was  judged  to  be  unacceptable  because  of  the  absence  of 
a  suitable  basement  area  (item  II),  and  although  the  specialty  stores  had  suitable 
basements,  the  building  did  not  have  sufficient  height  (item  I). 

6.  Lincoln  1st  Federal  Savings  Office  Building,  Spokane,  WA  (NSS  No.  248)  — 

This  structure  is  an  eight-story  building  constructed  in  1963.  The  first  two  levels 
are  reinforced  concrete  one-way  slab  and  beam  construction,  and  the  remaining 
floors  are  steel  frame  with  steel  beams  and  purlins  supporting  metal  deck,  topped 
with  concrete.  The  walls  are  of  metal  studs,  insulation,  and  various  interior 
finishes,  and  the  exterior  is  faced  with  marble.  The  building  has  two  full  levels 
below  grade. 

This  building  was  rejected  for  analysis  because  of  incomplete  drawings  -  only 
the  architectural  drawings  were  furnished  (item  III). 

Buildings  Scheduled  for  Demolition 

Another  source  for  selecting  candidate  buildings  was  obtained  through  our  work 
in  connection  with  FEMA  Contract  EMW-C-0582,  Building  Demolition.  It  was 
anticipated  that  a  minimum  of  eight  to  ten  buildings  would  be  explosively  demolished 
during  the  contract  period. 

The  use  of  these  structures  would  have  several  advantages.  They  are  more 
likely  to  be  completely  documented  structurally,  as  the  demolition  contractor  is 
required  to  make  a  thorough  investigation  in  order  to  determine  the  optimum  location 


for  placement  of  explosive  charges  (e.g.,  test  shots  are  normally  conducted  to 
determine  the  strength  of  component  parts)  and  they  are  usually  vacant  —  thus 
accessible  for  surveys.  Moreover,  the  post-demolition  data  would  be  valuable  as 
input  and  as  a  check  on  the  theoretical  analysis.  As  there  were  no  NSS  buildings 
suitable  from  the  plans  at  hand,  we  added  a  criterion;  i.e.,  item  VI  -  demolition  data 
should  be  available  to  compare  with  analysis. 

The  six  candidate  structures  that  would  possibly  be  available  included  the 
following: 


1.  Brewery  Bunding,  NJ  —  A  six-story  reinforced  concrete  structure.  This 
building  would  be  interesting  from  the  standpoint  of  an  experiment  in  that  the 
planning  was  to  demolish  only  the  top  floors  and  leave  the  lower  two  floors  intact  for 
future  use  as  a  shopping  mall. 

2.  Medical  School  Dormitory,  Minneapolis,  MN  —  A  seven-story  building  with 
floors  of  reinforced  concrete  pan  and  joist  construction  and  hollow  tile  walls.  The 
building  was  constructed  in  the  early  1930's. 

3.  Bank  Building,  Jacksonville,  FL  —  A  fifteen-story  building  constructed  of 
a  steel  frame  encased  in  tile. 

4.  Of  flee  Buildings,  Buffalo,  NY  —  Two  fifteen-story  buildings  constructed 
of  structural  steel. 

5.  Cornhusker  Hotel,  Lincoln,  NB  —  A  ten-story  hotel  constructed  in  1927  of 
reinforced  concrete  one-way  slab,  beam  and  girder  construction.  The  infill  walls 
were  hollow  clay  tile  covered  with  brick  veneer,  and  the  building  had  a  partial 
basement. 


6.  Olympic  National  Life  Building,  Seattle,  WA  —  A  twelve-story  reinforced 
concrete  3lab,  beam,  girder  and  column  building  constructed  in  1906.  The  walls  were 
reinforced  concrete  clad  with  sandstone  veneer.  The  building  had  a  full  basement. 
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Because  of  the  economic  climate,  during  the  past  year  only  Nos.  5  and  6,  the 
Cornhusker  Hotel  in  Lincoln  and  the  Olympic  National  Life  Building  in  Seattle,  have 
been  demolished  this  year.  Both  of  these  structures  were  very  old,  not  typical  of 
contemporary  construction,  and  it  was  therefore  felt  that  they  would  not  produce 
data  consistent  with  the  intent  of  the  contract.  Buildings  Nos.  1  through  4  were 
rejected  because  of  item  VI  —  there  would  be  no  demolition  data. 

Buildings  Previously  Demolished 

Further  sources  for  candidate  buildings  in  Contract  EMW-C-0582  included 
previously  demolished  buildings.  A  data  search  of  files  on  these  buildings  was 
performed  to  determine  if  any  buildings  previously  demolished  could  be  used.  Four 
candidate  structures  were  selected  for  analysis  from  the  files  of  Controlled 
Demolition,  Inc.  (CDI). 

1.  American  Industrial  Building,  Hartford,  CT  —  A  16-story  box  column 
steel-gusseted  structure.  This  was  an  interesting  structure  from  an  analysis  stand¬ 
point;  however,  no  plans  were  available  showing  the  structural  details  (item  III). 

2.  Abe  Lincoln  Hotel,  Springfield,  IL  —  A  14-story  structure,  cast-in-place 
concrete,  with  a  large  lobby  area.  Massive  concrete  beams  spanned  the  lobby,  and 
this  would  have  provided  an  excellent  example  of  the  "soft-story"  collapse 
mechanism;  however,  again  no  detailed  plans  were  available  (item  III). 

3.  Tutwiler  Hotel,  Birmingham,  AL  —  A  14-story,  steel  framed  structure. 
This  building  was  built  using  steel  produced  in  a  local  steel  mill,  and  was  designed 
with  the  exterior  columns  being  approximately  three  times  as  stiff  as  the  interior 
columns.  Thus,  this  building  was  eliminated  because  the  exterior  frame  was  not 
typical  of  most  buildings  (item  V). 

4.  Continental  Life  Building,  Atlanta,  GA  —  A  15-story,  cast-in-place 
reinforced  concrete  structure.  Constructed  about  1950,  the  building  covered  and 
area  90  feet  by  150  feet.  The  structure  had  a  full  basement  11  feet  high,  a  first 
floor  16  feet  high,  with  the  remaining  stories  being  10  feet  high  (see  Figures  2-2  and 
2-3). 


13 


Fig.  2-3.  Elevation  View  of  Continental  Life  Building,  Atlanta,  Georgia 


On  three  sides  the  wall  panels  were  of  masonry  and  glass,  while  the  fourth  side 
(on  the  north)  was  constructed  with  a  5-story  high  reinforced  concrete  shear  panel 
with  block  masonry  the  remainder  of  its  height. 

Supporting  columns  were  22  in.  square,  with  three  interior  rows  the  length  of 
the  building  supporting  one-way  pan  joist  beam  and  girder  construction.  Pans  were 
12  in.  by  20  in.  with  a  beam  length  of  22  ft  6  in.  and  a  12  in.  by  24  in.  structural 
section.  Girders  spanned  21  ft  6  in.  with  a  14  in.  by  28  in.  structural  section. 
Column  reinforcement  was  spiral  wrapped. 


BUILDING  SELECTION 

This  last  building,  the  Continental  Life  Building  (or  "Peachtree",  after  its 
street  address)  in  Atlanta,  met  nearly  all  of  the  criteria,  and  was  selected  for  the 
"program  tuning"  structure  under  this  contract.  The  required  data  on  the  as-built 
structure  itself,  as  well  as  information  on  its  demolition,  were  readily  available  from 
the  files  of  CDI,  the  demolition  contractor.  The  fact  that  this  building  was  a 
reinforced  concrete  structure,  and  thereby  focused  the  investigation  on  this  type  of 
structures,  does  not  in  any  way  restrict  the  generality  of  the  theoretical  analysis  or 
of  the  computer  modeling,  but  only  limits  the  breadth  of  the  verification. 
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Section  3 

TASK  2  -  UPGRADING  PLANS 


BASEMENT  UPGRADING  PLANS 

Using  the  guidance  developed  in  previous  SSI  reports  (Refs.  1  and  2),  and 
augmented  with  other  reports,  upgrading  plans  were  developed  for  the  structure 
selected.  Upgrading  plan  details  include: 

1.  Layout  of  the  basic  upgrading  system. 

2.  Selection  of  type,  size,  and  location  of  the  particular  structural 
members  used  in  the  upgrading.  Consideration  was  given  to  how  the 
material  would  be  brought  into  the  structure  and  installed. 

3.  Design  of  closures  to  seal  openings. 

4.  Special  upgrading  for  unusual  areas  in  or  near  the  designated  shelter 
space,  including  closures. 

5.  Specification  of  initial  radiation  and  fallout  protection. 

6.  Design  of  emergency  access  and  egress  structures. 

Experience  in  the  recent  MILL  RACE  test  event  (Ref.  27)  was  used  to 
supplement  items  1  through  4  in  the  upgrading  scheme. 

The  basement  area  for  upgrading  is  shown  in  plan  view  on  Figure  3-la,  and  the 
structural  details  of  the  first  floor  slab  are  given  in  Figures  3-lb  and  3-lc. 
Structural  upgrading  of  the  first  floor  is  necessary  to  enable  both  the  placement  of 
soil  to  provide  the  necessary  protection  against  radiation  and  to  provide  support  to 
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Plan  of  Basement  Showing  Column,  Elevator,  Stairway,  and 
Utility  Vault  Locations. 
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Fig.  3-lb.  Section  of  First  Floor  Showing  the  Framing  Plan. 
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the  floor  for  blast  overpressures  on  the  order  of  30  to  50  psi.  The  upgrading  plans, 
which  follow,  are  designed  for  a  protection  factor  (Pf)  of  1,000  against  fallout 
radiation. 

Radiation  Protection 

The  overhead  mass  required  for  radiation  protection  in  this  basement  shelter  is 
300  psf  to  achieve  the  Pf  of  1,000.  The  first  floor  slab  is  3  inches  thick  and 
provides  only  13%  of  the  required  overhead  mass.  Soil,  having  a  density  of  100  lb 
per  cubic  foot,  and  placed  to  a  depth  of  32  inches  over  the  entire  first  floor,  is  the 
most  expedient  solution  and  will  be  used  for  this  structural  upgrading.  The  basement 
area  90  feet  wide  by  150  feet  long  with  32  inches  of  soil  will  require  1,330  cubic 
yards  of  soil  to  be  hauled  to  the  site  and  placed  on  the  floor. 

Blast  Protection 

The  first  floor  must  be  upgraded  to  withstand  from  30  to  50  psi  blast 
overpressure  in  addition  to  the  weight  of  the  32  inches  of  soil,  which  is  required  for 
radiation  protection.  The  recommended  method  (Ref.  1)  for  a  floor  system  of  this 
type  is  to  use  post  and  beam  shoring  beneath  the  joists  and  post  shoring  beneath  the 
girders.  The  shoring  can  be  placed  no  farther  apart  than  the  quarter  span  distance 
for  the  members  being  shored.  Figures  3-2a  and  3-2b  show  the  post  and  beam 
shoring  that  is  required  beneath  the  joists,  and  Figure  3-3  shows  the  post  shoring 
that  is  used  beneath  the  girders. 

Shoring 

A  plan  of  the  post  and  beam  shoring  arrangement  for  the  basement  is  shown  in 
Figure  3-4.  It  should  be  noted  that  the  posts  and  beams  weigh  150  lb  and  350  lb, 
respectively;  the  timber  posts  and  beams  can  be  hauled  via  standard  elevators  into 
the  basement.  While  two  men  can  handle  and  erect  the  post  shoring,  the  post  and 
beam  shoring  will  require  mechanical  assistance.  Post  and  beam  shoring  can  be 
placed  by  two  men  with  the  aid  of  either  a  sheetrock  jack  or  a  forklift,  both  of 
which  are  available  at  most  equipment  rental  shops.  Actual  field  testing  has  been 
conducted  using  this  type  of  shoring  showing  its  feasibility  (see  Ref.  27).  In  order 
to  achieve  lateral  stability  of  the  entire  shoring  system,  it  is  recommended  that  the 
bracing  shown  in  Figure  3-5  be  installed. 


Fig.  3-2a.  Post  and  Beam  Shoring  of  Concrete  Joist 
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CLOSURE  DESIGN 


Two  c  'ferent  closures  will  be  required  for  the  basement:  the  first  involves 
the  sealing  off  of  the  two  stairways  leading  into  the  basement,  and  the  second 
requires  the  sealing  off  of  the  two  elevator  shafts.  The  stairs  can  be  sealed  off  by 
shoring  the  first  floor  landing  and  stairs  leading  into  the  basement  and  then  placing 
soil  to  a  depth  of  36  inches  above  the  first  floor  level.  Figure  3-6  shows  a  detail  of 
the  stairway,  shored  and  with  soil  placed  for  radiation  protection.  The  two  elevator 
shafts  prior  to  upgrading  should  have  the  elevators  run  up  to  about  mid  building 
height  so  that  both  the  elevator  cars,  counterweight,  and  associated  connecting 
cables  are  well  clear  of  the  basement  shelter  area.  Upgrading  of  the  elevator  shaft 
consists  of  placing  soil  to  a  depth  of  36  inches  above  the  first  floor  slab  (see  Figure 
3-7). 


Access  and  Ventilation 

Large  diameter  culvert  pipes  or  precast  concrete  manholes  could  be  placed  in 
the  stairway  or  elevator  shafts  prior  to  their  being  sealed  off,  and  having  soil  placed 
in  them  for  radiation  protection.  The  culvert  pipes  or  manholes  could  serve  the  dual 
purpose  of  providing  both  access  and  ventilation  for  the  shelter.  Additional 
ventilation  would  probably  be  obtainable  by  disconnecting  conduits  and  pipes  that 
connect  the  utility  vault  along  the  east  wall  of  the  basement. 
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Fig.  3-7.  Upgraded  Elevator  Shaft. 
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Section  4 

TASK  3  -  THEORETICAL  ANALYSIS 


INTRODUCTION 

This  task  is  the  heart  of  the  contract  work  and  treats  some  shelter  upgrading 
problems  for  the  first  time.  Some  work  has  been  performed  on  the  analysis  of  frame 
response  to  blast  and  earthquake  forces  and  the  failure  analysis  of  the  various 
components  such  as  columns,  beams,  and  floor  slabs  that  make  up  a  basement  shelter 
(see  Ref.  28);  but  little,  other  than  the  recent  SSI  work  to  develop  a  key  worker 
shelter  manual  (Refs.  1  and  2),  has  been  done  with  regard  to  the  analysis  of  frame 
response  on  an  upgraded  shelter.  Recent  full-scale  tests  of  a  number  of  upgraded 
reinforced  concrete  floor  systems  and  steel  frame  with  concrete  steel  deck  floor 
systems  indicate  that  the  modes  of  failure  for  upgraded  shelters  may  be  significantly 
different  from  those  for  non- upgraded  systems. 

It  was  found  that  the  classical  analysis  of  these  floor  systems,  based  on  simple 
or  continuous  span  ultimate  strength  design,  is  no  longer  a  valid  approach  in  failure 
prediction  when  considering  a  shored  system.  As  an  example  of  the  different  failure 
modes,  Fig.  4-1 A  shows  a  prestressed  precast  concrete  slab,  simply  supported  at  the 
ends,  which  was  tested  without  shores.  The  mode  of  failure  was  by  flexure,  as 
predicted.  The  slab  developed  positive  moment  tension  cracks  under  the  load  points 
early  in  the  test,  and  failure  occurred  in  flexure  at  these  locations.  This  type  of 
failure  was  predicted  by  conventional  design  methodology.  Fig.  4-1B  shows  one  end 
of  a  similar  slab  that  was  shored  at  midspan  and  loaded  to  failure.  The  failure 
shown  in  the  figure  occurred  as  a  result  of  bond  failure  of  the  prestressing  strands, 
causing  a  sudden  shear/flexure  near  the  end  support.  This  type  of  failure  was  not 
predicted  prior  to  the  test.  As  a  result  of  an  extensive  test  program  and  companion 
analytical  work,  SSI  has  developed  a  prediction  methodology  (Refs.  5  and  6)  by  which 
modes  of  failure  can  be  predicted  in  terms  of  the  locations  of  the  shores  and  the 
shear  and  moment  stresses  developed. 
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The  other  major  area  not  adequately  covered  by  past  research  is  the 
representation  of  the  complete  time  history  of  frame  collapse  behavior  under  blast 
loading;  this  includes:  the  elastic  phase,  the  elasto-plastic  formation  of  joint  yield 
hinges,  the  collapse  mechanism,  and  finally  the  fracture  of  elements  due  to  excess 
distortion.  The  dynamic  analysis  portions  of  the  STRUDL  program  and  the  concepts 
of  plastic  analysis  are  used  to  fulfill  this  research  need. 

The  complete  theoretical  effort  in  Task  3  was  divided  into  three  subtasks: 
3A  -  Development  of  loading  criteria 
3B  -  Development  of  joint  resistance  functions 
3C  -  Development  of  failure  modes 

SUBTASK  3A:  LOADING  FUNCTIONS 

Introduction 

The  blast  loading  on  a  structure  is  a  complex  function  of  the  incident  blast 
wave  (described  mainly  by  its  peak  overpressure  and  dynamic  pressure  wave  forms) 
interacting  with  structural  parameters  of  size,  shape,  orientation,  and  response.  To 
reduce  the  complexity  of  the  analysis .  of  the  interaction  process,  the  common 
approach  has  been  to  consider  a  structure  to  be  in  one  of  two  possible  categories, 
that  is,  either  as  a  diffraction-type  structure  or  a  drag-type  structure.  In  the 
former,  the  critical  response  is  to  the  peak  blast  wave  overpressure  while  in  the 
latter  the  response  is  to  the  entire  dynamic  pressure  pulse.  The  diffraction-type 
structure  of  principal  interest  here  is  a  large  building  with  strong  exterior  walls, 
typically  a  multistory  reinforced  concrete  building  with  small  window  area,  while 
typical  drag-type  structures  are  electrical,  radio  and  television  transmission  towers, 
and  truss  bridges. 

It  should  be  noted,  however,  that  few,  if  any,  true  diffraction-type  high  or 
medium  rise  structures  exist.  Research  in  Refs.  9,  10,  and  11  on  failure 
overpressures  for  various  wall  types  shows  that  most  nonreinforced  walls  fail  at  very 
low  overpressures,  less  than  5  psi.  Even  arched  nonreinforced  (not  normally  found  in 
high-rise  structures)  and  reinforced  concrete  walls  fail  at  overpressures  less  than  10 


to  12  psi.  Thus,  large  buildings  with  these  types  of  walls,  or  with  large  window 
areas,  have  been  commonly  considered  drag  structures  based  on  the  rationale  that 
brittle  wall  structures  can  resist  the  overpressure  for  only  a  fraction  of  the  duration 
of  the  blast  wave,  leaving  only  the  exposed  frame  to  bear  the  load.  While  it  is  true 
some  buildings  will  respond  mostly  as  drag-type  and  others  as  diffraction-type 
structures,  oversimplification  of  this  approach  ignores  the  real  situation  for  a  large 
number  of  buildings  wherein  frangible  walls  can  survive  long  enough  to  impart  very 
high  loads  into  supporting  frame  members  before  being  swept  out.  In  other  words, 
real  buildings  will  experience  both  types  of  loading  with  the  relative  importance  and 
contribution  of  each  type  in  causing  damage  depending  upon  wall  construction  and 
blast  wave  characteristics. 

An  indication  of  the  magnitude  of  the  load  that  can  be  imparted  to  the  frame 
by  frangible  portions  of  a  structure  was  obtained  from  some  shock  tunnel  tests  on  a 
wall  having  an  opening  both  with  and  without  glass  window  panes.  Data  from  two 
pairs  of  these  tests  are  shown  in  Figure  4-2.  It  will  be  noted  that  for  the  tests  at 
incident  pressure  of  p  =  3  psi,  the  load  on  the  frame  with  an  open  window  was 
approximately  50  kips,  but  with  window  glass  it  was  approximately  130  kips. 
Corresponding  differences,  140  kips  vs  250  kips,  were  noted  for  the  p  =  5.5  psi  tests. 

Figure  4-3  is  a  plot  of  the  impulse  transmitted  to  the  frame  by  an  unfailing 

wall  with  and  without  glass  in  a  27%  open  section.  The  net  difference  due  to  glass 

2  2 

is  (0.175  -  0.095)  lb-sec/in.  ,  or  0.080  lb-sec/in.  .  The  data  available  in  reduced 
form  are  limited  because  the  some  300  tests  done  in  the  shock  tunnel  have  never 
been  analyzed  for  impulse  (however,  some  data  were  presented  in  Ref.  29,  and  these 
are  shown  in  Figures  4-2  and  4-3);  there  was  no  interest,  until  now,  and  all  of  the 
shock  tunnel  data  are  at  much  lower  overpressures  than  are  of  interest  in  this 
program.  Therefore,  until  these  data  are  analyzed  and  some  higher  overpressure 
tests  conducted  in  a  field  test  (or,  perhaps,  scaled  in  a  small-scale  shock  tube)  it 
remains  necessary  at  present  to  make  a  best  estimate  of  loads  and  durations  of  the 
wall/frame  interactions  for  the  analytical  purposes  of  this  program.  A  discussion  of 
debris  and  loading  implications  is  presented  in  Appendix  C. 
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2.  Load  Data  from  Tests  With  and  Without  Glass  (27%  window  area). 
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Fig.  4-3.  Net  Impulse  Data  at  P  “12.7  psi  (P  ■  5.5  psi)  with  and 
Without  Glass  (27%  window  area) . 
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Development  of  Loading  Function 

The  technical  details  of  the  blast  loading  of  structures  are  taken  basically  from 
the  Effects  of  Nuclear  Weapons  (Ref.  30).  The  loading  condition  of  interest  is  that 
from  a  Mt  weapon  at  40  psi. 

For  a  closed  structure  the  front  face  loading  is  as  shown  below  where: 
p  =  peak  overpressure 

p(t)  =  time  variation  of  overpressure 

q  =  peak  dynamic  pressure 

q(t)  =  time  variation  of  dynamic  pressure 

pr  =  reflected  overpressure 

t+  =  duration  of  positive  phase  of  the  blast  wave 
ps  =  stagnation  pressure  =  p(t)  +  q(t) 
tg  =  clearing  time 


o  -fc,  TIME-  ■£* 

AVEUBASE  FRC2NJT  V+CJc.  UJADlNO 


When  the  blast  wave  strikes  a  flat  surface,  such  as  the  side  of  a  building  at  90 

degrees  (normal  incidence),  a  reflected  overpressure  occurs  which  is  more  than  twice 

the  incident  overpressure.  As  the  wave  front  passes  the  front  of  the  structure, 

rarefaction  waves  move  from  the  edges  across  the  front  of  the  structure  weakening 

the  loading  until  at  a  time,  t  ,  the  loading  has  been  reduced  to  the  stagnation  value. 

3  + 

The  stagnation  pressure  then  gradually  reduces  to  zero  at  a  time,  t  ,  as  the  incident 
loading  pulse  decreases.  The  time,  tg,  is  given  by  3S/U  where 
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S  =  the  height  or  one-half  width  of  the  structure 
U  =  the  shock  front  velocity. 


Values  of  pp,  q,  and  U  as  a  function  of  p  are  given  in  Figure  4-4  (Ref.  30). 

Once  the  shock  wave  has  engulfed  the  entire  structure, the  back  face  will  also 
be  loaded  as  illustrated  in  the  following  sketch.  This  assumes  that  the  pulse 
duration  is  very  long  compared  with  the  travel  time  down  the  structure  and  the 
clearing  times. 


The  net  horizontal  loading  is  that  given  by  subtracting  the  back  face  loading 
from  the  front  face  loading.  The  loading  pulse  shown  below  was  selected  as  a  first 
approximation  to  this  rather  complicated  pattern. 
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REFLECTED  PRESSURE  FOR  NORMAL  INCIDENCE 


PEAK  OVERPRESSURE  (PSD 


Fig.  4-4.  Relation  of  Ideal  Blast  Wave  Characteristics  at  the  Shock 
Front  to  Peak  Overpressure. 

Source:  Effects  of  Nuclear  Weapons,  Ref.  30. 


DYNAMIC  PRESSURE  (PSI ) 


For  a  structure  with  a  large  percent  of  its  surface  in  glass,  such  as  is  of  concern 
here,  Ref.  30  suggests  computing  the  loading  on  the  individual  structural  elements 
treating  them  as  closed  structures  and  then  summing  them.  This  does  not  change  the 
peak  loading,  which  is  still  the  peak  reflected  overpressure,  but  it  greatly  shortens 
the  duration  of  the  diffraction  spike,  since  the  half-width  of  the  structural  members 
is  used  for  S  in  place  of  the  half-width  of  the  entire  structure. 


From  the  sketch  of  a  typical  panel  of  the  structure  shown  below,  the  average  S 

value  is  1.1  ft,  so  that  t  is  1.6  msec. 

s 


This  time,  however,  completely  neglects  the  effects  of  the  windows,  which  as  pointed 
out  earlier  is  not  correct.  From  the  window  data  it  is  possible  to  calculate  an 
effective  clearing  time  as  follows:  the  impulse  transmitted  to  the  frame  from  the 
window  with  27%  glass  was  0.08  psi-sec;  extrapolating  this  to  the  68%  glass  case  of 
the  structure  of  concern  gives 

I  =  0.08  x  68/27  =  0.20  psi-sec 

Distributing  this  impulse  into  a  triangular  loading  spike  of  peak  value  equal  to  the 
peak  reflected  pressure  of  147  psi  gives  an  effective  clearing  time  of: 

t  =  21/p  -  (2  x  0.2)/l47  -  0.0027  sec 
w 

Adding  this  to  the  frame  clearing  time  of  1.6  msec  gives  a  total  of  4.3  msec,  which 
was  rounded  off  to  4  msec  for  the  computer  calculations. 
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For  the  drag  portion  of  the  loading  the  following  equation  was  used  for  the 
change  of  dynamic  pressure  with  time: 


t.\  ..  +v2  — 2(t/t  ) 

q(t)  =  q(l  -  t/t  )  e 


where  t  is  the  positive  phase  duration  of  overpressure.  For  the  time  range  of 
interest  ( 200  msec)  the  difference  in  impulse  between  this  equation  and  that  given 
by  the  data  in  Figure  4-5  (Ref.  30)  is  less  than  5%. 


A  drag  coefficient  of  2.0  was  used  in  calculating  the  loading  during  the  drag 
phase.  Ref.  30  suggests  a  value  of  1  with  the  loading  from  the  structural  elements 
of  both  the  front  and  back  faces  included.  In  the  present  calculation  only  the  front 
face  is  considered  so  that  a  large  drag  coefficient  is  appropriate. 


The  loading  characteristics  used  in  the  computer  calculations  are  summarized 
below  and  illustrated  in  Figure  4-6. 


Summary  of  Loading  Conditions 


Peak  overpressure 

40 

psi 

Weapon  yield 

Mt 

Positive  phase  duration  of  overpressure* 

2 

sec 

Peak  reflect  overpressure 

147 

psi 

Peak  shock  front  velocity 

2000 

ft/sec 

Peak  dynamic  pressure 

28 

psi 

Total  clearing  time 

4 

msec 

Total  area  of  panel 

215 

ft* 

Frame  area  of  panel 

69 

ft* 

Drag  coefficient 

2 

*  at  40  psi 
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NORMALIZED  DYNAMIC  PRESSURE,  q(f)/q 


Fig,  4-5.  Rate  of  Decay  of  Dynamic  Pressure  With  Time  for  Several 
Values  of  the  Dynamic  Pressure. 

Source:  Effects  of  Nuclear  Weapons,  Ref.  30. 
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It  should  be  noted  that  the  actual  shape  of  the  diffraction  loading  spike  as  well 
as  the  early  time  drag  loading  is  not  critical  since  the  response  of  the  structure  to 
this  loading  is  governed  by  the  total  impulse  delivered.  This  is  because  the  natural 
period  of  the  structure  (5  sec)  is  several  orders  of  magnitude  larger  than  the  loading 
duration.  Several  computer  analyses  were  performed  to  confirm  that  the  diffraction 
loading  is  indeed  an  impulse.  Figure  4-7  is  a  plot  of  mid-height  velocity  and 
displacement  resulting  from  the  4  msec  diffraction  loading.  This  impulse  alone  is 
probably  sufficient  to  collapse  the  structure,  as  the  displacement  at  0.5  sec  is 
sufficient  to  yield  the  lower  beams,  and  the  velocity  greater  than  over  12  in.  per 
second. 


SUBTASK  3B:  JOINT  RESISTANCE  FUNCTIONS 

Introduction 

The  analysis  of  building  structures  under  blast  loading  requires  a  knowledge  of 
the  effective  elastic  limit  capacity,  M^,  for  the  beam  and  column  joints.  This 
section  provides  methods  of  evaluating  these  capacities  for  the  structural  elements  in 
both  steel  and  reinforced  concrete  frame  buildings.  The  developments  are  limited  to 
the  beam  and  column  elements  of  frames,  since  the  common  brittle  exterior  and 
interior  wall  construction  will  be  removed  by  the  peak  overpressure  of  the  30  to  50 
psi  incident  blast  wave.  Further,  one  principal  objective  of  the  analysis  is  to 
determine  the  effect  of  the  frame  column  failure  on  the  ground  floor  (shelter  ceiling) 
slab.  The  presence  of  structural  walls  would  complicate  any  findings  concerning  the 
integrity  or  failure  of  this  slab  due  to  the  continuous  column  interaction.  Both  the 
steel  and  the  reinforced  concrete  member  capacities  will  be  evaluated  at  the 
strength  (ultimate  strength  design)  basis. 

Structural  Steel  Construction 

The  AISC  Manual  (Ref.  31)  provides  strength  information  for  beams,  columns, 
welded  and  bolted  connections,  and  splices.  Part  2  of  this  manual  gives  the  specific 
plastic  design  or  strength  values  (See  interaction  diagram,  Figure  4-8,  and  Table  4-1 
for  all  definitions): 


44 


TABLE  4-1:  AISC  MANUAL,  PART  2  DEFINITIONS 


m 


P 

P. 


Gross  area  of  an  axially  loaded  compression  member 

Coefficient  applied  to  bending  term  in  interaction  formula  for  prismatic  members  and 
dependent  upon  column  curvature  caused  by  applied  moments 


Fa  Axial  compressive  stress  permitted  in  a  prismatic  member  in  the  absence  of  bending 
moment  (kips  per  square  inch) 

Fy  Specified  minimum  yield  stress  of  the  type  of  steel  being  used  (kips  per  square  inch) 
M  Factored  bending  moment  (kip-feet) 

M  Critical  moment  that  can  be  resisted  by  a  plastically  designed  member  in  the  absence  of 
axial  load  (kip-feet) 


Plastic  moment  (kip-feet) 
Factored  axial  load  (kips) 
Euler  buckling  load  (kips) 


cr  Maximum  strength  of  an  axially  loaded  compression  member  or  beam  (kips) 
S  Section  modulus 
Z  Plastic  section  modulus 
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Beams, 


M  =  ZF  , 

p  y . 

or  Mp  =  SFy  for  older  non-compact  sections. 


Columns, 

P/P  +  (CM)/  f(l  -  P/P  )  M  "1  <  1 
cr  m  L  e  m  j  — 

where  P  =  1.7  AF  ,  C  =  0.4  for  reversed  curvature,  M  =  M  . 
cr  a  m  m  p 

(Strengths  for  connections  and  splices  are  1.7  times  the  corresponding  AISC  Manual 
Part  1  allowable  stress  values.) 


With  respect  to  beam  flexural  capacities,  it  is  quite  probable  that  older 
construction  may  have  steel  (builtup  or  rolled)  sections  that  are  non-compact,  such 
that  they  would  buckle  before  developing  full  plastic  capacity.  For  these  sections 
the  M  value  should  be  taken  as  SF  ,  for  the  tabulated  or  calculated  section  modulus, 


In  any  given  structure,  and  particularly  in  older  structures,  it  may  be  possible 

that  column  splice  details  (such  as  shown  in  Figure  4-9)  may  constitute  a  weak  link 

in  column  M  values.  If  these  splices  are  weak  in  flexural  resistance  and  near  the 
P 

column  base,  then  their  estimated  Mp  capacity  should  be  used  at  these  locations. 
The  transverse  shear  resistance  of  the  splice  may  also  be  a  weak  link  and  should  be 
investigated. 


Also,  particularly  in  older  construction,  the  interior  beam  column  connections 
may  be  simple  web  or  flange  clip  angles  with  rather  minimal  M^  values  (see  Figure  4- 
10).  This  type  of  detail  should  be  identified,  and  the  appropriate  estimated  value 
should  be  used  in  the  analysis. 


Recommended  ductility  ratios,  y  =  M  /M  are  8  to  10  for  fully  developed 

r  p 

sections  or  rigid  connections,  where  Mp  is  the  failure  moment.  Other  than  for  the 
case  of  weak  splice  details,  the  shear  capacity  of  a  steel  section  will  not  be  less  than 
the  shear  necessary  to  develop  the  M  value.. 
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Moment  diagram 


Moment  diagram 


Reinforced  Concrete  Frame  Construction 

By  far  the  most  prevalent  forms  of  construction  are  the  flat  slab  (or  waffle 
slab)  and  the  two-way  slab  systems,  see  Figure  4- 11 A  and  B.  The  CRSI  Handbook 
(Ref.  32)  provides  most  of  the  strength  information  required  for  these  systems  along 
with  the  beam  and  column  section  capacities  for  frame  elements. 


For  the  majority  of  slab  system  frame  structures,  the  equivalent  frame  is 
defined  by  a  frame  strip  along  each  column  line,  having  a  width  equal  to  the  bay 
width  perpendicular  to  the  column  line.  Figure  4-12  shows  the  general  flexural 
section  configurations.  The  flexural  strength  of  these  beam  sections  can  be 
taken  from  the  CRSI  Handbook  using  the  assumption  that  strengths  for  F^  =  60  ksi 
steel  with  the  <f>  factor  are  equal  to  strengths  for  F  =  40  ksi  without  the  <f>  factor. 
The  most  realistic  estimate  of  would  be  without  the  <f>  factor  multiplier. 


For  a  suitable  approximate  estimate  of  beam  capacities,  M  ,  the  following 

P 

procedure  can  be  applied  to  the  construction  elements  in  Figure  4-12. 


o  Assumption:  Positive  steel  area  equals  one-third  negative  steel  area  at  the 
column  face 

Negative  M  =  (Negative  A  )F  (0.8h) 
u  s  y 

Negative  A  =  top  steel 
s 

Positive  M  =  1/3  Negative  M 
u  B  u 

o  Assumption:  Negative  dead  load  moment  equals  one-third  Negative  (see 
Figure  4-13) 

Negative  M  =  Negative  M  -  M__  =  2/3  M 
p  °  u  DL  u 

Positive  M  =  1/3  Negative  M  +  M_,  =  2/3  M 
p  &  u  DL  u 

Therefore,  Negative  M  =  Positive  M 

P  P 
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A.  The  Flat  Slab 


B.  The  Two-Way  Slab 


Fig.  4-11.  Examples  of  Slab  Systems. 


51 


NEGATIVE.  Mu  <£APA£ITY 


APPROXIMATE 


Rtsitive  Mu  CAPja^rry 
=  '/3>  Ne^Arive.  Mu 


Mp  REi-^ricPM^S. 


Fig.  4-13.  Approximate  M  Relations. 


For  column  M  values,  under  the  effects  of  axial  and  flexural  load  interaction, 
P 

the  interaction  tables  of  the  CRSI  Handbook  provide  the  M  =  P  times  the 

P  u 

eccentricity  e.  These  values  are  for  F  =  60  ksi  with  the  4>  factor  and  can  be 
assumed  to  be  equal  to  the  for.  F^  =  40  ksi  without  the  <t>  factor.  Figure  4-14 
shows  an  example  interaction  curve,  and  Figure  4-15  shows  the  CRSI  table  values  for 
the  Peachtree  Building  columns. 


It  is  important  to  recognize  that  columns  may  fail  in  shear;  this  capacity  can  be 
estimated  by 

V  =  A  250  A  lb 

u  c  c  c 

where  A  is  column  section  area  in  square  inches. 

c  j- 

The  shear  stress  of  4'T  represents  the  presence  of  shear  stirrup  steel  in  the  form  of 
column  ties. 


Also,  reinforcing  steel  splices  may  provide  a  weak  link  if  the  spliced  bars  are 
not  staggered;  one-half  might  be  used  when  this  splice  condition  is  present. 


Recommended  ductility  ratio  w  =  M  /M  is  5  for  fully  developed  sections  or 

F  p 

rigid  frame  joints. 


Joints  at  the  Ground  Level  Slab  Intersection 

Although  the  joints  that  exist  at  the  ground  level;  i.e.,  at  the  superstructure  of 
the  basement,  are  similar  to  other  joints  in  the  structure,  they  behave  much 
differently.  The  reason(s)  for  this  behavioral  difference  is  that  the  upgrading  of  the 
basement  structure  will  greatly  change  the  floor  stiffness.  The  slab  portion  of  the 
structure  will  be  shored,  which  will  increase  its  stiffness  relative  to  the  as-built, 
pre-upgrading  stiffness.  The  dead  load  will  also  be  increased  two  to  three  times 
owing  to  the  depth  of  soil  needed  for  radiation  protection.  In  addition,  the  slab 
restricts  lateral  motion  at  the  ground  level  if  the  basement  is  a  satisfactory  shelter 
space,  i.e.,  the  slab  and  basement  walls  must  be  poured  integrally.  Typical 
upgrading  schemes  that  illustrate  these  differences  are  shown  in  Figure  4-16. 
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Fig.  4-15.  Column  Interaction  Table  for  Peachtree  Building. 
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Close-in  Shoring 

For  cases  where  there  are  high  moments  induced  in  the  ground-level  slab; 
either  due  to  heavy-debris  loads  or  due  to  the  "crow-bar"  prying  action  of  strong 
steel  frame  columns,  the  shores  may  be  placed  "close-in"  to  the  columns  in  order  to 
assist  the  slab  or  a  weak  steel  frame  connection  in  resisting  the  imposed  column 
movement,  see  Figure  4-17. 


SUBTASK  3C:  DEVELOPMENT  OF  FAILURE  MODES 

The  Role  of  Plastic  Analysis  Techniques 

In  the  introductory  paragraph  of  this  section,  the  comment  about  collapsing 
frames  constituting  a  "new  frontier"  turned  out  to  be  truer  than  even  the  authors 
initially  perceived.  In  the  late  50's  and  early  60's,  the  plastic  design  of  framed 
structures  was  a  primary  part  of  structural  engineering.  The  concept  of  plastic 
design  had  the  potential  for  economies  in  building  design  and  construction  and  for  the 
development  of  a  more  uniform  factor  of  safety  throughout  the  structure  (i.e.,  a 
more  uniform  probability  of  collapse  of  the  various  components).  However,  the 
evolution  of  plastic  design  of  framed  structures  has  almost  ceased  since  the  recent 
advent  of  sophisticated  computer  programs.  These  computer  programs  permitted 
rapid  economical  linear  elastic  analysis  of  even  the  most  complex  frame  structures, 
and  hence,  the  simple  "plastic  mechanism"  analysis  advantage  of  plastic  design  has 
been  outmoded.  Also,  the  material  economies  of  plastic  design,  inherent  in  member 
section  capacities  and  in  moment  redistribution,  have  been  incorporated  into  building 
code  provisions.  Hence,  present  design  is  performed  by  use  of  elastic  analysis  for 
stresses,  and  section  design  by  approximations  of  plastic  capacity. 

For  the  purposes  of  the  frame  response  history  analysis,  however,  there  was  a 
specific  need  to  identify  the  different  mechanisms  of  frame  collapse  and  the 
respective  collapse  loads.  This  "mechanism"  analysis  is  based  on  the  resistances  of 
the  beams  and  columns  for  a  typical  frame  as  determined  by  the  methods  in  Subtask 
3B  and  on  the  particular  computer  model.  The  form  of  analysis  and  prediction  of 
failure  modes  in  this  section  follows  the  classical  plastic  design  techniques  presented 
in  Appendix  A.  By  the  use  of  the  computer,  the  fundamental  theories  of  limit  design 
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Fig.  4-17.  Shoring  Countermeasure  Against  High  Induced  Moments  in 
Ground-Level  Slab. 
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will  predict  the  most  probable  failure  modes  for  a  building  frame.  The  initial  fully 
elastic  response  of  the  framed  structure,  however,  may  show  some  anomalies,  in  that 
higher  moments  may  be  achieved  at  levels  other  than  those  predicted  by  the  limit 
design  failure  mode  analysis.  These  moments  will  be  elastic  failures,  however,  and 
as  successive  triads  approach  the  mechanism  of  collapse,  classical  plastic  design  will 
eventually  dominate  the  failure  mode.  This  will  be  demonstrated  in  Section  5, 
Computer  Analysis. 


Figure  4-18  shows  a  representation  of  a  typical  framed  structure  model.  The 

hinge  formation  in  this  type  of  structure,  if  the  basement  is  properly  shored,  will  be 

forced  in  the  column  above  the  floor  level.  At  the  column-beam  intersection  at 

joint  "A",  in  the  comer  at  the  first  level,  as  long  as  the  column  is  continuous  through 

the  joint  and  the  beam  has  some  moment  reistance,  there  should  be  yield  hinge 

development  in  the  column.  Because  of  the  relative  stiffness  of  the  floor  system,  it 

is  very  possible  that  in  the  elastic  analysis  initial  yielding  may  develop  in  the  beam 

at  prior  to  complete  hinge  formation  in  the  column  labeled  M  .  However,  since 

side  sway  is  prohibited  at  this  ground  level,  a  hinge  will  ultimately  develop  at  M  , 

and  only  small  damage  will  occur  in  the  area.  For  the  interior  joint  at  nB",  ^ 

occurs  twice,  M  is  on  the  basement  side  of  the  column,  and  M  is  on  the  column 
c  p 

above  the  first  story.  Here  again,  if  the  column  is  continuous  through  the  joint,  and 
has  any  value  at  all,  the  hinge  will  form  at  In  concrete  frames  this  would 
nearly  always  be  the  case,  since  the  column  steel  is  generally  spliced  above  the  floor 
level,  thereby  creating  a  zone  of  weakness. 


As  discussed  in  Subtask  3B,  many  steel  frames  have  essentially  pin  connections 
or  seated  connections  at  the  beam-to-column  connection,  and  the  column  is  indeed 
continuous  through  the  joint.  This  results  in  a  lack  of  moment  capacity  in  the  beam, 
but  this  deficiency  may  be  overcome  by  locating  shores  near  the  columns.  Placing  of 
shores  near  a  column  generates  significant  moment  resistance  in  a  member,  even  if  it 
is  nominally  pin  connected.  This  resistance  occurs  because  the  shear  capacity  of  the 
joint  times  the  short  lever  arm  distance  to  the  shore  generates  sufficient  moment  to 
cause  the  hinge  to  occur  above  the  floor  level. 
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For  reasons  discussed  previously,  the  basement  selected  should  include  a 
reinforced  concrete  floor  cast  monolithically  with  the  wall  system  in  order  to  resist 
side  sway  at  ground  level. 

Simplified  Plastic  Analysis  of.  Peachtree  Building 

The  particular  frame  used  as  our  example,  and  which  will  be  discussed  below  in 

more  detail,  is  a  small  portion  out  of  the  larger  Peachtree  Building,  the  structure 

that  is  analyzed  in  Section  5.  The  moment  capacities  for  the  beams  and  columns  in 

this  frame  are  based  on  the  calculations  shown  in  Subtask  3B.  The  moment  capacity 

of  the  beam,  M  ,  is  approximately  200  kip-ft  for  both  positive  and  negative  moment 
P 

directions,  and  the  moment  capacity  of  the  columns  in  the  lower  story  is 
approximately  483  kip-ft,  or  2.41  of  the  beams. 

Figure  4-19  is  a  sketch  of  the  basic  frame  used  in  this  mechanism  study.  It  is 
a  three-story,  two-bay  frame  with  a  basement  structure  restrained  at  the  ground 
level.  The  frame  consists  of  15  joints,  10  reactions  external,  and  20  bars,  and 
results  in  15  independent  mechanisms  of  collapse.  One  interesting  observation  is 
that  all  of  the  possible  beam  mechanisms  (see  Appendix  A)  are  not  developed  as  a 
result  of  side  loading.  This  differs  from  what  is  usually  presented  in  the  typical 
plastic  design  text,  where  often  gravity  loads  are  the  prime  concern  in  building 
design.  In  this  study,  the  blast  load  dominates,  shown  as  P  on  this  structure,  and  the 
gravity  loads  are  merely  the  dead  loads  of  the  structure.  Hence,  in  the  limit 
analysis  for  these  structures,  only  the  lateral  loads  are  considered;  the  dead  load  has 
been  nominally  taken  care  of  in  the  adjustment  of  the  moment  capacities  for  the 
various  beam  and  column  components. 

Another  very  interesting  aspect  of  this  particular  problem  is  that  of  the  15 
independent  mechanisms,  12  are  degenerate  modes  or  joint  rotation  modes  where 
external  work  is  not  done,  and  this  leaves  only  three  independent  side-sway  modes 
where  external  work  is  done.  This  implies  that  far  more  complex  problems  can  be 
performed  by  hand  analysis  than  was  originally  thought  possible.  In  a  sense,  there 
are  only  three  independent  mechanisms,  and  these  can  be  used  to  very  closely 
approximate  the  true  failure  load. 
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Fig.  4-19.  Side-Sway,  Hinge,  and  Combined  Mode  Assessment. 
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Fig.  4-19.  (contd) 
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Fig.  4-19.  (contd) 
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The  three  independent  side-sway  modes,  as  illustrated  in  Figure  4-19,  are  first 
level  (1),  second  level  (2),  and  third  or  top  level  (3),  respectively.  Combined  Mode 
1  is  the  combination  of  side-sway  Mode  1,  side-sway  Mode  2,  plus  the  appropriate 
independent  hinge  modes.  Combined  Mode  2  is  a  combination  of  all  three  side-sway 
modes,  and  again,  the  appropriate  independent  hinge  modes.  Although  other  possible 
modes  of  failure  exist,  it  is  felt  that  these  are  the  dominant  modes  based  on  the 
simplicity  of  the  blast  loading  system. 

Table  4-2  is  a  summary  of  the  analysis  of  the  frame  described  above  with  beam 
strength  to  column  strength  ratios  representing  strong  columns,  medium  columns,  and 
weak  columns.  The  first  case,  strong  columns,  could  occur  where  design  criteria 
require  substantial  resistance  to  lateral  loads,  such  as  in  earthquake  zones,  and  the 
third  case,  weak  columns,  is  not  uncommon  when  the  span  between  the  columns  is 
fairly  long  compared  with  the  story  height,  or  the  vertical  loads  are  very  heavy. 
The  medium  column  case  takes  into  account  a  structure  lying  between  these  two 
extremes. 


TABLE  4-2:  SUMMARY  OF  FRAME  ANALYSIS 


M  ,  /H  Ratio 
column  beam 

(1)  (2)  (3) 


MODE  Strong  Column  Medium  Column  Weak  Column 

2.41  M  /M  »T~M  /M  M  /vT“m 

_ £ _ E _ 2 _ E _ P  P 


Independent 

1 

3.62 

Mp/2 

*2.12 

Mp/2 

*1.50  Mp/2 

Independent 

2 

9.64 

v* 

5.66 

Mp/2 

4.00  Mp/2 

Independent 

3 

22.5 

Mp/2 

15.30 

Mp/2 

12.00  Mp/2 

Combined  1 

3.34 

Mp/2 

2.26 

Mp/2 

2.11  Mp/2 

Combined  2 

*3.30 

Mp/2 

2.61 

Mp/2 

3.02  Mp/2 

*  indicates  minimum  or  failure  P. 
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As  shown  in  Table  4-2,  the  Combined  Mode  2  is  the  collapse,  or  failure,  mode 
when  columns  are  considerably  stronger  than  the  beams.  In  this  case,  the  entire 
building  leans  over,  and  failures  occur  at  the  column  bases,  just  above  the  floor 
grade,  and  then  lateral  movement  of  the  upper  stories  occurs  similar  to  a  falling  tree. 
In  the  case  of  columns  that  are  only  slightly  stronger  than  beams,  the  shift  is  toward 
Independent  Mode  1,  where  the  upper  stories  remain  intact,  and  the  building  moves 
sideways  above  the  first  level.  It  should  be  noted  that  the  differences  between 
loads  for  Independent  Mode  1  and  Combined  Mode  2  are  not  great.  Accordingly, 
where  structures  have  a  fairly  equal  beam/column  resistance,  the  modes  can  be  either 
"soft  story"  with  the  upper  stories  remaining  rigid,  or  "falling  tree". 

In  the  third  analysis,  where  the  beam  is  stronger  than  the  columns,  there  is  a 
very  large  difference  between  the  Independent  Mode  1  and  Combined  Mode  2.  As 
the  structure  fails  in  Mode  1,  it  moves  laterally  at  the  first  story  and  basically 
remains  a  vertical  and  intact  system  above.  This  type  of  "soft  story"  is  often  noted 
in  earthquake  catastrophes,  such  as  the  failure  of  Olive  View  Hospital  in  San 
Fernando  and  several  of  the  apartment  buildings  in  Caracas,  Venezuela.  In  these 
cases,  the  structures  had  been  designed  as  rigid  frames,  but  above  the  first  story, 
masonry  infill  walls  had  actually  formed  a  shear  wall  building  such  that  the  frame  did 
not  respond  in  a  frame-like  manner,  but  as  a  rigid  box  on  a  one-story  frame.  It  is 
important  to  note  that  a  similar  phenomenon  would  occur  in  a  true  frame  by  simply 
changing  the  column-to-beam  strength  ratio. 

Another  observation  that  can  be  made  from  this  analysis  is  that,  in  many  steel 
frame  buildings,  the  exterior  frames  are  designed  to  carry  all  the  lateral  load,  and 
the  interior  frames  are  "pinned"  at  the  beam- to- column  connection;  i.e.,  the  interior 
frames  are  actually  in  the  so-called  combination  Mode  3  by  design.  Hence,  one  may 
expect  these  steel  frames  to  be  forced  into  Combined  Mode  2  even  though  the 
moment -resisting  frame  in  the  building  itself  may  not  otherwise  be  governed  by  this 
mode.  As  an  example,  if  the  Peachtree  Building  had  been  a  steel  frame  building,  the 
interior  five  frames  may  have  been  framed  in  the  manner  of  Combined  Mode  2,  with 
all  the  beams  pinned  to  the  columns,  and  the  two  end  frames  would  have  probably 
been  designed  as  rigid  frames  with  continuous  beams  and  columns  in  order  to  be  the 
moment  resisting  or  lateral  load  resisting  system  of  the  whole  structure.  In  modern 
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building  construction,  the  floors  are  generally  poured  concrete  on  steel  decking  and 
act  as  a  diaphragm  to  carry  the  load  to  the  rigid  frames.  When  there  are  several 
frames  in  the  center  of  the  building  that  are  actually  constructed  in  this  manner, 
Combined  failure  Mode  2,  there  is  a  high  probability  that  the  moment  resisting  frames 
at  the  end  would  be  forced  into  this  mode  shape,  and  the  calculated  lower  bound 
failure  mode  may  be  exceeded. 

A  simplified  analysis  of  the  Peachtree  Building  is  summarized  in  Figure  4-20. 

It  examines  the  entire  15  stories  and  4  bays  in  brief.  Because  it  is  known  that  in 

this  particular  structure,  the  column  strength  is  much  greater  than  the  beam 

strength,  one  would  expect  the  combined  modes  could  be  important.  Figure  4-20 

illustrates  the  various  modes  of  failure  and  provides  the  corresponding  value  of  the 

load,  P,  from  the  simplified  analysis.  It  is  apparent  from  the  latter  that  the  lower 

bound  static  failure  mode  is  Combined  Mode  5,  where  P  =  0.73  M  /i  (a  static  load  of 

P 

14.6  kips).  This  represents  hinges  in  four  floors  of  beams  and  hinges  in  the  columns 
at  the  base  and  below  the  sixth  floor.  The  static-elastic  computer  analysis  predicts 
first  hinges  forming  at  P  =  10  kips  in  the  second  story  beams,  and  if  beams  did  not 
fail,  columns  would  fail  at  about  P  =  17  kips,  which  of  course  bounds  our  plastic 
lower  bound  of  14.6  kips. 

It  is  interesting  to  note  that  the  predictions  for  the  various  failure  modes 
(represented  by  the  P  values  in  the  figure)  are  relatively  close;  i.e.,  ±15%  from  the 
average  modal  value,  which  is  really  quite  close,  considering  our  knowledge  of  the 
blast  phenomena  and  actual  construction  details  as  built.  Further,  one  would  expect 
the  response  under  dynamic  loads  to  migrate  toward  lower  modes  (perhaps  as  low  as 
Mode  2  or  Mode  3)  once  the  contribution  of  the  upper  story  interiors  is  taken  into 
account.  This  added  resistance  to  motion  is  due  to  inertial  forces  of  the  upper 
stories,  which  tend  to  make  the  upper  stories  lag  (see  Section  5,  Figure  5-10). 

Another  consideration  would  be  reduction  of  column  strength,  as  one  goes  up  in 
the  building.  The  effect  here  is  that,  if  the  column  bending  strength  is  reduced  as 
one  moves  to  upper  stories  (because  of  the  reduced  vertical  load  requirements  in  the 
building),  this  would  reduce  the  variation  in  P  values  above  the  fourth  floor  and 
hence,  the  collapse  mode  in  one  of  the  higher  side-sway  modes.  Even  so,  it  is 
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Fig.  4-20.  (contd) 


unlikely  that  this  effect  could  cause  a  collapse  at  a  significantly  lower  load  with 
typical  building  design. 

For  example,  in  the  Peachtree  Building  it  would  be  difficult  to  reduce  the 
lateral  resistance  (or  the  bending  resistance)  of  the  columns  below  about  50%  base 
column  and  maintain  compliance  with  building  code  requirements.  The  beam 
strengths  would  be  expected  to  remain  constant  throughout  the  structure,  so  that  in 
the  upper  stories  one  could  visualize  a  change  in  tendency  toward  a  combined  mode 
effect  in  the  upper  regions  and  the  side-sway  mode  in  the  lower  regions  of  the 
building.  An  assessment  of  the  above  conditions  shows  such  a  change  would  reduce 

Mode  14  from  0.96  M  /*•  to  about  0.93  M  /«.. 

P  P 

Variations  such  as  these  should  be  checked  in  future  work,  and  some  simple 
bounding  rules  need  to  be  evolved  for  quick  gross  analyses  of  high-rise  buildings. 
Another  important  observation  is  that,  for  this  large  scale  building,  the  difference 
between  the  side-sway  mode  of  the  first  story  and  the  Combined  Mode  2  is  only  16%; 
although  the  columns  are  much  stronger  than  the  beams  in  this  particular  building, 
the  condition  for  a  column-induced  failure  and  a  beam-induced  failure  is  relatively 
small.  Hence,  a  possible  approach  to  use  in  order  to  bound  failure  loads  for  the 
shelter  space  is  to  look  at  the  two  most  divergent  modes  and  assume: 

(1)  that  the  building  collapse  is  in  a  mode  like  the  side-sway  (soft  story) 
mode,  with  the  first  floor  columns  falling  and  the  building  moving  as  a  rigid  body 
laterally  until  the  upper  stories  collapse  on  the  basement  somewhat  as  a  unit; 

(2)  that  the  building  collapse  is  in  a  combined  mode,  where  the  beams  and 
columns  fail  and  all  or  part  of  the  building  rotates  like  a  falling  tree. 

Consequences  and  effects  on  the  shelter  space  will  be  investigated  for  these  two 
cases. 


From  the  results  of  the  simplified  analyses  (Figure  4-20)  the  constant  or 
statically  applied  (not  time  varying)  failure  load  P  can  be  determined,  or  at  least 
bracketed  to  an  accuracy  of  about  25%.  For  the  Peachtree  Building  this  failure 
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blast  load  is  P  =  0.87  M  /£,  op  about  17.4  kips  (vs  14.6  statically).  This  is  the  result 

P 

of  the  blast  pressure  at  each  story  level.  When  this  static  value  of  P  is  compared 
with  the  time- varying  blast  level  used  in  the  computer  analysis  of  Section  5,  it  is 
quite  obvious  that  failure  will  occur  at  a  very  early  time. 

Objectives  of  the  Dynamic  Analysis 

Knowing  that  the  "static"  failure  load,  P,  is  on  the  order  of  15  kips  and 
comparing  this  with  the  very  high  blast  load,  P  shown  in  Figure  4-21  (a  straight-line 
approximation  to  Figure  4-6),  it  is  quite  obvious  that  frame  failure  will  occur. 
Therefore,  the  main  questions  to  be  addressed  are: 

(1)  What  is  the  primary  mode  of  failure,  "soft  story"  or  "falling  tree"? 

(2)  Is  the  blast  load  intensity  and  duration  sufficient  to  cause  progressive 
mechanisms  in  the  structure  after  the  formation  of  the  primary  failure  mode? 

(3)  Will  the  entire  structure  be  virtually  shattered  before  it  drops  to  the 
ground,  or  will  it  remain  intact  and  destroy  itself  by  gravity  action  and  impact  on  the 
shelter  slab? 

(4)  Will  the  imparted  velocity  and  residual  blast  force  be  sufficient  to  carry 
the  failing  structure  beyond  the  shelter  area,  or  will  most  of  the  structure  come 
down  on  the  shelter  slab? 

In  Section  5,  the  complete  successive  elastic,  elasto-plastic,  plastic  mechanism 
analysis  is  performed  for  the  given  blast  loading  of  Subtask  3A  and  member  strength 
properties  of  Subtask  3B. 
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4-21.  Comparison  of  Dynamic  Blast  Load  and  "Static"  Failure  Load. 
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Section  5 

COMPUTER  MODEL  DEVELOPMENT, 
ANALYSIS,  AND  RESULTS 


INTRODUCTION 

It  has  been  common  practice  in  civil  engineering  to  represent  time-varying 
dynamic  loads  as  equivalent,  pseudo-static  loads  for  the  analysis  of  structures.  This 
applies  to  blast,  wind  and  earthquake  loading  conditions.  However,  as  more  data  are 
obtained  regarding  the  nature  and  effects  of  transient  dynamic  loadings,  structural 
engineers  are  becoming  increasingly  aware  of  the  necessity  to  perform  more  realistic 
dynamic  analysis  in  order  to  predict  structural  response. 

Also,  acceptable  procedures  for  generating  pseudo-static  loadings  have  not 
Jt>«6n  developed  for  many  dynamic  loading  conditions.  The  requirement  for  dynamic 
analyses  leads  to  a  direct  need  by  the  engineer  for  a  sophisticated  general  purpose 
computer  software  system.  For  the  analysis  in  this  report  the  static  and  dynamic 
features  of  the  Georgia  Tech  Integrated  Civil-  Engineering  System  Structural  Design 
Language  (GTSTRUDL)  were  selected. 

Figure  5-1  shows  the  dynamic  capabilities  of  the  chosen  computer  software 
system.  Several  types  of  dynamic  analysis  are  available.  The  program  offers  a 
choice  of  performing  a  shock  analysis  (response  spectrum)  to  determine  the  maximum 
effect  due  to  an  input  frequency  spectrum,  or  the  choice  of  a  normal  mode  analysis, 
or  a  direct  integration  of  the  equations  of  motion  to  determine  the  transient  response 
due  to  specified  transient  loads  or  support  accelerations.  A  description  of  the 
fundamental  computational  methods  in  the  dynamic  analysis  program  is  presented  in 
Appendix  B. 
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GTSTRUDL  Dynamics  Overview 


DESCRIPTION  OP  BUILDING  STRUCTURES  FOR  ANALYSES 


Two  separate  structures  were  analyzed  under  both  static  and  dynamic  blast 
loading  conditions.  A  three-story,  two-bay  structure  was  used  to  verify  the  method 
and  accuracy  of  the  analysis  procedure.  A  fifteen-story,  four-bay  concrete 
structure  with  basement*  was  used  as  the  representative  structure  to  predict  the 
collapse  mechanism  when  subjected  to  a  30  to  50  psi  blast  pressure  level.  To 
facilitate  the  analysis  and  reduce  the  amount  of  data  to  be  handled,  only  one 
representative  cross  section  through  the  full  width  and  height  of  the  building  was 
used.  Verification  of  the  chosen  computer  code  consisted  of  modeling  the  three- 
story  building  using  15  nodes,  12  column  elements,  and  8  beam  elements.  Total 
degrees  of  freedom  numbered  36,  considering  the  structure  as  a  planar-moment 
resisting  frame  ignoring  the  out-of-plane  displacements  and  rotations.  Boundary 
conditions  consisted  of  the  columns  being  fixed  at  ground  level  and  the  first  floor 
being  resisted  from  displacing  horizontally. 

Static  analyses  performed  on  the  structure  were  verified  using  moment 
distribution  techniques.  Dynamic  analyses  were  verified  using  approximations  to 
fundamental  frequencies  and  mode  shapes.  A  single  mass-spring  system  was  used  to 
verify  time  behavior  of  the  structure.  The  computer  code  applied  checked  out  well 
against  the  behavior  predicted  via  the  simplified  analysis.  However,  to  lend 
credence  to  the  collapse  predictions  for  the  representative  structure  actual  data  are 
needed. 

The  representative  structure  was  modeled  for  static  analysis  using  101  nodes, 
80  column  elements,  72  beam  elements  and  8  basement  shoring  elements.  The  column 
elements  were  22  in.  square  reinforced  concrete,  beam  elements  were  12  in.  x  24  in. 
reinforced  concrete,  and  basement  shoring  elements  were  3  in.  x  8  in.  wood  posts. 
(See  Figure  5-2  and  5-3  for  elevation  and  floor  plan.)  All  beam  and  column 
intersections  were  considered  moment  resisting,  except  that  the  wood  shoring  was 
considered  pin  connected  to  the  basement  floor  and  concrete  first  floor  beams.  In 
addition,  for  the  dynamic  analyses  a  lumped  mass  system  was  used  using  88  masses 
with  a  total  of  237  dynamic  degrees  of  freedom. 


# 


Peachtree  Building,  demolished  by  CDI  in  1980,  located  in  Atlanta,  Georgia. 


For  dynamic  history  analysis,  a  modal  analysis  method  of  solution  was  used,  as 
described  in  the  appendix  under  the  heading  of  computer  analysis  fundamentals. 
Figures  5-4,  5-5,  and  5-6  show  the  joint  and  member  numbering  scheme,  member 
forces  at  any  time  instant,  and  the  shape  of  the  deformed  structure  under  static  load. 


RESULTS  OF  ANALYSES  OF  PEACHTREE  BUILDING 

To  determine  the  behavior  of  the  structure  when  subjected  to  a  blast,  a 
conversion  must  be  found  to  translate  the  blast  into  a  force  history  as  it  impinges  on 
the  structure.  Subtask  3A  developed  this  conversion,  and  Figure  5-7  shows  the 
resulting  load  function  as  it  applies  to  the  joints  on  the  side  of  the  building  facing 
the  blast.  The  application  of  the  blast  load  in  addition  to  the  dead  load  of  the 
building  are  shown  pictorially  in  Figure  5-8.  Joint  resistances  were  developed  in 
Subtask  3C,  the  column  resistance  is  480  kip-feet  and  the  beam  resistance  200  kip- 
feet  to  reach  the  yield  limit  of  the  concrete. 

For  the  pseudo-plastic  phase,  the  reduced  stiffness  E  l  was  taken  as  one-tenth 
of  the  initial  elastic  stiffness  El.  The  ductility  factor  for  the  ratio  of  failure 
(fracture)  rotation  to  elastic  limit  rotation  D  was  taken  as  five.  Ordinary 
reinforced  concrete  construction  could  have  a  range  of  equal  to  3  to  6,  and  the 
selected  value  gives  the  model  a  chance  to  exhibit  a  reasonable  amount  of  mechanism 
development  before  failure. 

Initial  failure  due  to  yielding  under  blast  load  conditions  occurs  60  milliseconds 
(0.06  seconds)  after  initial  blast  impingement,  when  a  collapse  mechnism  as  shown  in 
Figure  5-9  has  formed.  Total  collapse  of  the  structure  occurs  when  all  the  yield 
hinges  have  reached  their  ultimate  moment  carrying  capacity.  Failure  of  the 
collapse  mechanism  starts  at  120  milliseconds  and  total  failure  is  complete  at  200 
milliseconds  after  initial  blast  impingement. 

The  lower  two  floors  have  completely  collapsed  and  are  in  the  process  of 
depositing  rubble  within  a  few  feet  of  their  original  position.  The  rest  of  the 
structure  is  free  from  its  foundation  and  is  now  a  free  falling  object  governed  by  its 


86 


J  iiiPPUlf  JOHI  ^QCAllSNS 


Fig.  5-4.  Member  and  Joint  Numbering  Scheme. 
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Fig.  5-5.  Example  Computer  Output  for  Member  Stresses. 
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Fig.  5-6.  Deformed  Structure  Static  Analysis. 
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5-7.  Blast  Load  Input  Function. 
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Fig.  5-9.  Collapse  Mechanism. 
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mass,  velocity,  acceleration  and  applied  forces.  Figures  5-10  and  5-11  show  the 
horizontal  displacement  of  the  structure  vs  time.  The  lower  two  floors  are  severely 
distorted  while  the  upper  floors  move  more  or  less  as  a  unit  without  appreciable 
structural  damage.  During  this  phase  the  only  damage  mechanism  to  the  shelter  is 
the  blast  loading  itself  as  the  columns  above  are  yielding.  Figure  5-12  shows  a 
graphic  history  of  each  phase  of  the  structure's  behavior. 


COMPUTER  TIME  HISTORY  FOR  SHELTER  ELEMENTS 

The  more  detailed  results  of  the  computer  analysis  provide  the  history  of  the 
individual  member  actions  in  a  typical  bay  of  the  basement  shelter.  Referring  back 
to  the  general  member  code  or  numbering  scheme  in  Figure  5-4,  the  load-time  curves 
for  the  columns,  shore  struts,  and  shelter  ceiling  beams  are  given  in  Figures  5-13,  5- 
14,  5-15,  and  5-18.  It  should  be  noted  that  these  member  reactions  do  not  include 
the  effect  of  the  40  psi  vertical  overpressure  loading  on  the  shelter  ceiling.  From 
the  previous  approximate  analysis,  this  40  psi  pressure  does  not  fail  the  slab.  It  is, 
however,  effective  in  canceling  out  the  indicated  tension  forces  in  the  shoring  struts. 
(These  were  modeled  as  pinned-end  links  for  the  analysis,  but  in  reality  they  could 
not  take  tension.)  The  computer  indicated  30  kips,  and  this  is  completely  canceled 
by  an  order  of  150  kips  compression  due  to  the  blast  overpressure. 


SUMMARY  OF  ANALYSIS  CONCERNING  EFFECTS  ON  SHELTER  SPACE 

There  are  two  important  results,  one  very  encouraging  with  respect  to  the 
safety  of  the  basement  shelter  area  and  the  other,  rather  discouraging.  Taking  the 
encouraging  aspect  first,  the  failure  mode  of  the  frame  is  by  yield  mechanisms  in  the 
first  and  second  story  and  then  fracture  or  collapse  of  the  first  story  columns  at 
their  ductility  limit.  The  failure  moments,  axial  loads,  and  shears  developed  by 
these  first  story  columns  we  not  sufficient  to  cause  any  significant  damage  to  the 
beams,  slabs,  and  columns  of  the  upgraded  basement  shelter  area.  All  of  these 
shelter  members,  and  the  shoring  members,  have  reserve  strength  at  the  point  of  the 
first  story  frame  failure. 
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Fig.  5-10.  Early  Time  History  (to  t  =  0.12  sec) 
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Fig.  5-11.  Extended  Time  History  (to  t  =*  0.5  sec) 
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Fig.  5-14.  Representative  Member  Time  Histories. 
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The  discouraging  aspect  of  shelter  survivability  is  as  follows:  the  first  and 
second  story  mechanisms,  and  the  first  story  column  failures,  take  place  very  rapidly 
without  time  for  the  upper  stories  to  settle  down  "softly"  onto  the  shelter  slab.  For 
the  very  probable  case  of  brittle  column  failure,  it  is  quite  possible  that  the  entire 
building  could  have  virtually  free  fall  for  at  least  a  height  equal  to  the  first  story. 
This  falling  building  mass  might  be  cushioned  somewhat  by  the  failure  mechanisms  of 
the  first  and  second  stories.  The  impact  zone  could  be  over  most  of  the  shelter  area 
since  the  blast  force  and  its  imparted  velocity  to  the  building  at  the  time  of  collapse 
is  not  sufficient  to  carry  the  building  beyond  the  shelter  plan  area.  The  impact 
force  of  the  nearly  intact  upper  story  portion  of  the  building  could  well  exceed  the 
upgraded  shelter  slab  capacity.  This  result  definitely  needs  further  verification,  for 
other  frames,  and  for  three-dimensional  frame  models,  and  from  field  observations 
obtainable  from  building  demolition  experience.  Some  representative  numerical 
calculations  are  given  below. 


DEBRIS  DEPTHS  AND  BASEMENT  AREA  SURVIVAL 

The  frame  response  depicted  in  Figure  5-12  indicates  the  formation  of  plastic 
hinges  in  approximately  60  milliseconds.  Total  time  frame  for  building  failure  and 
collapse  of  the  first  two  stories  occurs  at  0.2  seconds.  The  loading  impulse  on  the 
building  is  sustained  for  4  seconds,  causing  the  building  to  collapse  both  vertically 
and  horizontally. 

As  stated  in  the  preceding  section,  the  shelter  survives  the  lateral  displacement 
and  failure  of  the  columns  at  the  first  floor  plastic  hinge  locations  (see  Figure  5-17). 
In  addition,  the  shelter  basement  ceiling  is  impacted  by  the  falling  debris  from  the 
building  collapse  and  must  sustain  the  weight  of  the  debris  plus  the  related  impact 
load  forces.  The  15  stories  result  in  a  total  debris  weight  at  an  equivalent  120  psf 
for  each  floor  and  12  psf  for  contents: 

132  psf  x  15  stories  =  1,980  psf,  or  14  psi 

Impact  of  the  falling  debris  is  estimated  at  two  times  the  static  weight  of  the  debris, 
or  28  psi. 
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A  rough  estimate  of  punching  capacity  of  the  shored  floor  system  is  240  kips 
per  shore,  and  an  8x8  shore  would  not  fail  below  approximately  300  kips.  If  the 
impact  load  from  the  falling  debris  occurred  in  conjunction  with  the  total  static 
debris  load  (not  likely),  total  load  on  a  shore  would  be  approximately  180  kips,  and 
therefore  not  cause  punching  failure.  Note  that  this  loading  does  not  include  any 
punching  effects  that  could  be  due  to  falling  columns  impacting  vertically  on  the 
slab.  This  type  of  punching  has  not  been  observed  in  our  experience  with  the 
demolition  of  concrete  frame  buildings;  however,  it  does  occur  in  steel  frames. 

Total  collapse  time  for  a  15-story  building  is  on  the  order  of  7  to  10  seconds. 
Thus,  the  blast  wave  forces  acting  on  the  shelter  ceiling  have  passed  prior  to  impact 
of  the  building  debris  or  the  debris  static  weight. 

The  recovery  of  the  shelterees  is  another  matter.  Based  on  our  experience  in 
building  demolition  with  CDI,  expected  debris  depths  are  20  inches  (0.5  m)  per  story. 
Thus,  15  stories'  debris  depth  would  approximate  an  equivalent  7.5  meters,  or  24.5 
feet.  This  debris  would  result  from  both  the  building  analyzed  and  from  all  adjacent 
buildings  of  equivalent  height. 

Also,  the  environmental  problems  of  dust  and  other  air  contamination  effects 
after  the  blast  event  are  not  covered  by  this  report  and  definitely  need  further 
study. 
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Section  6 

SUMMARY  AND  CONCLUSIONS 


This  report  presents  the  results  of  a  program  to  develop  a  theoretical  analysis  of 
the  effects  of  frame  response  on  basement  shelters  in  tall  buildings.  The  objective 
was  to  determine  the  effect  on  a  upgraded  basement  key  worker  shelter  of  the 
aboveground  portion  of  the  structure  being  subjected  to  a  blast  wave  (30  to  50  psi) 
that  would  destroy  the  building. 

This  program  investigated  both  steel  and  reinforced  concrete  frame  structures 
with  the  most  emphasis  on  poured-in-place  reinforced  concrete  beam,  slab,  and 
girder  type  framing  and  the  poured-in-place  flat-slab  and  flat-plate  type  of 
construction.  These  types  are  very  common  in  the  National  Shelter  Survey  inventory 
of  upgradable  structures. 

A  prediction  technique  was  developed  using  both  hand  and  computer  analysis. 
This  technique  was  tested  using  a  previously  explosively  demolished  15-story  cast-in¬ 
place  reinforced  structure,  the  Continental  Life  Building  in  Atlanta,  Georgia.  The 
results  of  this  analysis  indicated  that  the  upgraded  basement  would  have  survived 
even  though  the  aboveground  portion  of  the  structure  was  exposed  to  50  psi. 

Large  portions  of  the  debris  landed  on  the  shelter  roof  so  that  serious  questions 
remain  as  to  the  advisability  of  using  such  structures  as  shelters  because  of  the 
possible  problems  of  entrapment  of  the  shelterees.  Other  problems  that  need  to  be 
addressed  are  the  large  amounts  of  dust  created  in  the  collapse  process  and  the 
possibility  of  fire  in  the  debris  pile. 

It  also  should  be  noted  that  the  primary  work  in  this  program  was  on  reinforced 
concrete  structures,  where  the  punching  effects  of  falling  columns  do  not  appear  to 
be  a  problem.  Based  on  the  experiences  of  an  explosive  demolition  contractor, 
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Controlled  Demolition,  Inc.,  this  appears  to  be  a  serious  problem  in  steel  framed 
structures.  It  is  recommended  that  future  work  in  this  area  thoroughly  investigate 
the  collapse  of  steel  frames  both  theoretically  and  in  conjunction  with  building 
demolitions.  Also,  as  was  noted  in  Section  4,  very  little  information  is  available  on 
the  load  imparted  to  a  frame  by  failing  walls  at  pressures  above  15  psi.  Tests  at  the 
higher  pressures  of  interest  (30  to  50  psi)  need  to  be  conducted. 
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APPENDIX  A 

APPLICABLE  CONCEPTS  OF  PLASTIC  ANALYSIS 


APPLICABLE  CONCEPTS  OF  PLASTIC  ANALYSIS 


INTRODUCTION 

In  order  to  best  show  how  plastic  analysis  can  be  applied  to  frame  response,  it 
is  appropriate  to  introduce  some  of  the  basic  concepts.  Ref.  A-l,  "The  Plastic 
Design  of  Steel  Frames"  by  Lynn  S.  Beedle,  published  in  1958,  is  one  of  the 
pioneering  works  on  plastic  design  in  framed  structures  and  is  still  one  of  the 
fundamental  books  on  this  subject.  This  text  will  be  used  as  one  of  the  principal 
references  in  this  discussion.  The  primary  concept  used  in  plastic  design  is  taking 
advantage  of  the  ductility  of  a  material.  Figure  A-l  is  an  idealized  plot  of  stress  vs 
strain  for  steel.  When  elastic  design  concepts  are  used  with  a  typical  steel 
structure  the  allowable  stresses  are  on  the  order  of  20,000  psi  for  the  A7  steel; 
obviously  when  a  structure  is  operating  under  service  load  conditions  it  is  desirable 
to  keep  the  stress  at  or  below  such  a  stress  level.  However,  when  a  structure  is 
loaded  up  to  the  yield  stress,  34.1  ksi  in  this  particular  case,  it  does  not  mean  that 
the  structure  is  going  to  collapse,  but  merely  that  the  steel  has  reached  its  yield 
stress;  one  may  continue  to  stretch  or  distort  the  structure  and  it  will  still  remain 
intact,  but  will  be  permanently  distorted.  The  area  under  the  stress-strain  curve  is 
a  measure  of  the  energy  absorption  capacity  of  the  material;  as  can  be  seen  from  this 
figure,  the  energy  under  the  elastic  portion  is  small  compared  to  the  area  under  the 
entire  curve.  In  fact  the  curve  shown  in  Figure  A-l  is  only  the  first  10%  of  the 
entire  stress-strain  curve.  The  entire  curve  is  plotted  on  Figure  A-2  and  puts  some 
perspective  on  the  tremendous  ductility  of  steel. 

By  way  of  contrast,  Figure  A-3  shows  the  stress-strain  curve  for  plain  concrete 
in  compression.  Here  an  ultimate  compression  strain  of  about  0.3%  occurs  as 
opposed  to  a  value  of  20%  or  more  for  steel.  It  should  be  mentioned  that,  although 
concrete  is  far  less  ductile  than  steel,  in  real  world  structures  the  full  ductility  of 
steel  is  generally  not  realized.  The  joints  are  assembled  with  bolts  and  welds,  there 
are  dimensional  restraints,  and  the  gross  distortions  are  generally  limited  to  the 
region  of  the  curve  shown  in  Figure  A-l.  That  is,  perhaps  the  first  1%  to  3%  of  the 
distortion  is  really  available  for  energy  absorption  of  a  structure  prior  to  failure. 
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Steel  and  reinforced  concrete  beam  behavior  is  shown  by  the  curves  in  Figure 
A-4.  Here  M  is  the  applied  moment,  M  is  the  yield  moment  for  that  particular  cross 
section  or  beam,  is  the  ultimate  plastic  moment,  4>  is  the  rotation  angle  of  the 
joint  in  the  zone  of  plastic  deformation  in  the  beam  and  <t>  is  rotation  at  yielding  of 
the  outer  fiber  of  a  steel  beam  or  of  the  reinforcing  steel  in  the  concrete  beam. 
The  particular  plot  shows  two  curves  for  steel,  one  for  a  round  cross  section  and  one 
for  a  wide-flange  cross  section,  and  the  ratio  of  M  to  M  .  This  is  known  as  the 

P  y 

shape  factor;  typically  in  rectangular  sections  the  shape  factor  is  1.5,  and  in  wide- 
flange  sections  it  is  equal  to  1.12  to  1.14,  or  a  12%  to  14%  increase  from  yield 
moment  to  ultimate  plastic  moment  in  capacity  in  a  wide-flange  section.  The  <|>  to 
axis  is  a  measure  of  ductility  and  in  a  typical  reinforced  concrete  beam  the 
ductility  is  usually  in  the  range  of  3  to  6,  sometimes  greater  for  doubly  reinforced, 
carefully  designed  beams;  but  in  most  buildings  encountered  across  the  country  this 
range  of  3  to  6  for  ductility  is  not  typical.  The  effect  of  this  plastic  or  nonlinear 
section  behavior  on  the  total  structure  behavior  will  be  discussed  below. 


The  following  definitions  are  important  to  the  understanding  of  plastic  analysis: 

(1)  A  plastic  hinge  is  a  zone  of  yielding  due  to  flexure  in  a  structural 
member.  Although  its  length  depends  on  the  geometry  and 
loading,  in  most  of  the  analytical  work  it  is  assumed  that  all 
plastic  rotation  occurs  at  a  point.  At  those  sections  where 
plastic  hinges  are  located,  the  member  acts  as  if  it  were  hinged, 

except  with  a  constant  restraining  moment,  M  . 

P 

(2)  Plastic  hinges  form  at  points  of  maximum  moment.  Thus  in  a 
framed  structure  with  prismatic  members,  it  would  be  possible  for 
plastic  hinges  to  form  at  point  of  concentrated  load,  at  the  end  of 
each  member  meeting  at  a  connection  involving  a  change  in 
geometry,  and  at  the  point  of  zero  shear  in  a  span  under 
distributed  load. 


(3)  The  plastic  moment  equals  F^Z,  where  Z  is  the  plastic  section 
modulus. 


(4)  The  shape  factor  (f  =  Z/S)  is  one  measure  of  reserve  strength 


beyond  the  elastic  limit  (S  =  section  modulus). 
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Fig.  A- 
Source : 
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.  Typical  Concrete  Stress-Strain  Curve,  Short-Term  Loading. 
McCormac,  "Design  of  Reinforced  Concrete,"  Ref.  A-2. 


.  The  Maximum  Strength  of  a  Round  Pin  Compared  With  That  of 
a  Wide  Flange  Beam. 

,  "Plastic  Design  of  Steel  Frames,"  Ref.  A-l. 


Source:  Beedle 


The  effect  of  plastic  hinging  on  the  behavior  of  structures  can  most  easily  be 
illustrated  by  a  simple  example,  see  Figure  A-5.  As  the  load  P  is  applied  to  this 
simple  beam  at  midspan,  the  load  would  increase  in  a  linear  manner  until  the  outer 
fiber  reaches  yield.  This  would  be  in  the  range  of  classical  elastic  design  procedure. 
As  the  structure  is  continually  loaded  beyond  the  yield  point,  it  does  not  fail  and 
collapse,  but  the  yielding  progresses  from  the  outer  fibers  of  the  beam  section 
toward  the  centroid  of  the  beam  until  the  entire  section  has  either  gone  into  the 
plastic  domain  or  fracture  occurs,  as  in  the  case  of  reinforced  concrete.  The  upper 
picture  of  the  rectangular  section  in  Figure  A-5  can  also  be  represented  as  the 
normalized  M— <f>  curve,  Figure  A-6.  Figure  A-6  is  the  normalized  M/M^  vs  <j>/4> 
curve  for  a  rectangular  section  of  a  steel  beam.  This,  however,  is  very  nearly 
identical  to  that  of  reinforced  concrete.  The  point  labeled  l'  on  the  curve  for  the 
case  of  concrete  would  be  the  cracking  moment  in  concrete  when  the  first  initial 
tension  cracks  form;  point  (^2  on  the  curve  would  be  the  value  commonly  used  in 
the  design  of  concrete  structures  prior  to  applying  a  safety  factor;  and  point  V)  in  a 
concrete  beam  would  be  analogous  to  the  so-called  M  ,  or  the  ultimate  strength  of 
the  concrete  beam  near  =4.  Point  4J1  illustrates  the  immense  ductility  of  a 
steel  beam  in  this  mode  of  bending,  whereas  the  concrete  beam  will  more  than  likely 
fracture  somewhere  between  points  (3)  and  I\  Individual  resistance  functions  or 

M  values  are  described  under  Subtask  3B  (see  Section  4). 

P 

Figure  A-7  is  a  plot  of  M  versus  $  for  a  typical  wide-flange  section  of  steel. 
Here  note  that  the  shape  factor  is  small;  that  is,  there  is  only  a  14%  difference 
between  point  of  yielding  and  the  point  of  full  plastic  cross  section. 


The  full  impact  of  the  plastic  hinge  concept  is  not  really  appreciated  without 
the  idea  of  a  mechanism  of  failure  and/or  the  moment  redistribution  that  occurs  in  a 
continuous  or  indeterminate  structure  when  it  is  loaded  to  failure.  Figures  A-8  and 
A-9  illustrate  a  simple  beam  with  fixed  ends,  uniformly  loaded.  As  the  load  is 
increased,  the  first  fibers  to  yield  will  be  at  the  ends  of  the  beam  A  and  B;  these  are 
shown  in  Figure  A-8  as  ends  A-B  on  the  M-<t>  plots.  At  the  time  that  they  are  in 
this  idealized  picture,  they  provide  a  constant  plastic  moment,  but  the  center  of  the 
structure  is  still  well  into  the  elastic  domain;  that  is,  its  stresses  are  less  than  one 
half  the  way  to  yield,  hence  there  is  no  danger  whatsoever  of  the  structure 


A-5 


Fig.  A-5.  Theoretical  Length  of  Yielded  Portion  of  (a)  Rectangular 
and  (,b)  WF  Beam  With  Central  Concentrated  Load. 

Source:  Beedle,  "Plastic  Design  of  Steel  Frames,"  Ref.  A-l. 


Fig.  A-8.  Redistribution  of  Moment  in  a  Fixed -Ended  Beam  With 
Uniformly  Distributed  Load. 


Fig.  A-9.  Moment  Diagram  at  Various  Stages  of  Loading  for  Fixed-Ended 
Beam  With  Uniformly  Distributed  Load. 


Source:  Beedle,  "Plastic  Design  of  Steel  Frames,"  Ref.  A-l. 


collapsing.  It  will,  therefore,  require  a  continual  increase  of  load  to  develop  a  hinge 
in  the  middle  of  the  beam  in  order  to  achieve  collapse;  i.e.,  a  mechanism  shown  by 
the  heavy  dark  lines  of  (c)  in  Figure  A-8.  Figure  A-9  shows  the  sequential  increase 
of  moment  from  the  elastic  to  the  first  yield  at  the  end,  to  the  first  yield  at  the 
center,  to  the  development  of  a  mechanism  or  hinges  at  all  three  points  and  collapse. 
This  collapse  load  is  approximately  1.52  times  yield  load  for  a  wide  flange  section 
and  1.85  for  a  solid  rectangular  section. 

It  is  interesting  to  compare  these  results  with  elastic  theory  where  one  would 
predicate  failure  load  at  first  yielding,  when  in  fact  for  a  wide-flange  beam  to 
collapse  requires  52%  more  load  beyond  yield.  By  contrast,  if  a  rectangular  section 
were  bent  and  taken  to  collapse,  one  would  have  to  increase  the  load  85%  beyond 
yield  to  develop  full  plastic  hinges  at  all  three  points  and  thus  develop  a  mechanism 
of  failure  or  collapse.  Since  most  civil  structures  or  buildings  are  made  with  wide- 
flange  beams,  the  maximum  reserve  capacity  one  would  generally  encounter  in  a  beam 
of  this  type  would  be  the  52%  number.  In  fact  this  is  seldom  realized,  in  that  the 
beam  connections  have  limitations,  and  as  will  be  seen  in  more  complex  structures, 
there  is  a  smoothing  or  a  tendency  of  the  structure  to  behave  in  the  most  expedient 
way  to  collapse;  i.e.,  to  form  a  lower  bound  on  collapse  strength  nearer  to  that 
predicted  by  elastic  theory.  In  concrete  structures  the  extra  strength  from  first 
yield  to  a  collapse  mechanism  is  not  large.  This  is  because  the  concrete  moment 
capacity  is  controlled  by  steel,  so  that  the  steel  placed  in  a  concrete  beam  is 
selected  to  match  or  nearly  match  the  actual  load  moment  value  at  that  section; 
hence  there  is  little  excess  strength,  but  greater  redundancy  in  reinforced  concrete- 
structures.  In  both  concrete  and  steel  structures  the  code  attempts  to  reduce  the 
overdesign  by  allowing  a  redistribution  of  moments;  that  is,  the  code  allows  one  to 
allocate  a  portion  of  the  high  negative  moment  to  the  positive  moment  region,  thus 
reducing  excess  structural  capacity,  and  approximating  formation  of  a  mechanism. 

To  this  point  the  ideas  surrounding  plastic  behavior  of  structures  have  been 
discussed.  In  order  to  increase  the  potential  for  this  type  of  analysis  it  is  important 
to  add  some  theorems  of  plastic  analysis.  Formal  proof,  discussion,  etc.,  concerning 
these  theorems  can  be  found  in  several  references  such  as  Hodge  (Ref.  A-3)  and  Moy 
(Ref.  A-4). 
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The  Uniqueness  Theorem:  For  any  set  of  conditions  of  loads  and  strengths 
there  is  but  one  way  (one  mode)  to  collapse  the  structure  and  satisfy  the  four 
following  conditions. 

1)  Kinematic- Conditions:  Enough  plastic  hinges  to  form  a  mechanism  of 
collapse. 

2)  Equilibrium:  Statics  must  be  satisfied  for  all  parts  of  the  loaded 
structure. 

3)  Yield  Limits  must  be  observed;  the  absolute  value  of  the  moment  must  be 
equal  to  or  less  than  capacity  everywhere  in  the  structure  (M  _i  M^). 

4)  The  Dissipation  Condition:  There  must  be  positive  work  at  every  plastic 
hinge.  The  direction  of  the  load  bending  moment  must  be  in  the 
direction  of  rotation. 


Upper  Bound  Theorem:  is,  simply  stated,  the  violation  of  Item  3  of  the 
Uniqueness  Theorem,  i.e.,  exceeding  the  moment  M  >  M  (overshoot)  somewhere  in 
in  the  structure. 

Lower  Bound  Theorem:  is  the  violation  of  Item  1  of  the  Uniqueness  Theorem, 
that  is,  insufficient  number  of  hinges  to  form  the  kinematic  conditions  for  collapse. 

Independent  Mechanisms:  If  p  is  the  number  of  possible  plastic  hinges,  r  is  the 
redundancy  or  degree  of  indeterminancy  of  the  structure,  and  m  is  the  number  of 
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independent  mechanisms,  then 


m  =  p  -  r 


where  r  =  3n  -  3b  -  n  -s 

r 

and  n  =  number  of  joints 

b  =  number  of  members 

n  =  number  of  support  conditions,  and 
r 

s  =  special  conditions  such  as  built-in  hinges 
For  a  more  complete  discussion,  see  Wilbur  and  Norris  (Ref.  A-5). 

A  two-span  beam  with  plastic  section  capacity  will  be  used  to  illustrate  the 
use  of  these  theorems,  see  Figures  A-10  and  A-ll. 


Consider  Mode  1  (see  Figure  A-ll): 

Work  in  =  Work  out 

P(l)  =  M  (♦  +  +_)  +  M 

p  A  D  P 


Uniqueness  (1) 
Uniqueness  (2) 
Uniqueness  (3) 
Uniqueness  (4) 


kinematics  satisfied 
to  be  verified 
to  be  verified 
dissipation  satisfied 


Solving 


P  =  M  (2/£.  +  2/1)  +  M  2  H 
P  P 


P  =  6M/i 
P 
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MODE  1 


Moment  diagram  for  Mode  1 


MODE  2 

=  2/1,  ,j>  =  1/J l 

B  c 


Moment  diagram  for  Mode  2 


"POSSIBLE"  MODE  3 


Possible  moment  diagram 
for  Mode  3 


Fig.  A-ll.  Failure  Modes  for  Two-Soan  Beam  Example. 


This  is  either  a  solution  or  an  upper  bound. 


To  satisfy  equilibrium  the  moment  diagram  must  be  as  shown 


Hence  we  have  an  Upper  Bound  solution  and  have  violated  condition  (3),  i.e., 
exceeded  capacity. 


Upper  bound  solution 


The  lower  bound  solution  can  be  found  by  reducing  all  moments  by  3/4,  and  this 

creates  a  violation  of  condition  (1)  for  kinematics.  Hence,  we  know  then  the  true 

solution  is  greater  than 

3/4(6  M  /*)  =  9/2(M  / 1) 

P  P 

and  we  know 

9/2(M  h)  <  P  <  6M  fl 
.  P  P 

from  a  single  solution  of  the  problem. 


Consider  Mode  2  (see  Figure  A-ll): 

work  in  =  work  out 

P(l)  =  M  (  *  +  *  )  +  M  (<f>  _) 

p  □  C  p  o 

P  =  M  (2/2  +  l/£)  +  M  2  U 
P  P 


P  =  5M^/2,  which  is  either  an  upper  bound  or  the  solution. 


Equilibrium  Moment  Diagram 


Therefore  the  Mode  2  solution  satisfies 
Uniqueness  (1)  -  kinematics 
Uniqueness  (2)  -  equilibrium 
Uniqueness  (3)  -  |  M] _<  M^  everywhere 
Uniqueness  (4)  -  dissipation 

Therefore  P  =  5M  /«.  is  the  solution. 

P 

Consider  "possible"  Mode  3  (see  Figure  A-ll): 

The  uniqueness  theorem  part  (4)  requires  a  positive  work  or  energy  dissipation. 


+  <£  ,  +M  +M,  -  <p 

Hence  Mode  3  is  inadmissible,  as  it  violates  the  uniqueness  theorem. 

The  power  of  the  foregoing  theorems  and  arguments  for  our  particular  problem 
of  "Frame  Collapse"  is  very  great  indeed.  Any  assumed  mode  that  satisfies  (1) 
kinematics  and  (4)  dissipation,  will  provide  the  upper  and  lower  bound  solutions. 
Further,  if  these  bounds  are  reasonably  close  there  is  little  rational  reason  for 
attempting  to  find  the  "exact"  solution.  The  loading  conditions  under  blast  effects 
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have  a  wide  range  of  possibilities  and  variations.  The  actual  resistance  functions 
(plastic  moment  and  ductility)  have  great  variability,  plus  the  unknowns  of  actual 
construction  and  design  variabilities  are  also  present.  All  of  these  variabilities  lead 
to  the  realization  that  reasonably  close  bounds  can  provide  a  satisfactory  estimate  of 
the  collapse  load. 

To  illustrate  the  foregoing,  the  following  frame  is  presented. 


Hence 

h  =  10  hinges 

r  =  5  redundancies 

m  =  5  independent  mechanisms 

The  five  independent  mechanisms  are  shown  in  Figure  A- 12. 
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The  Solution  (see  Figure  A-12) 


(?)  q  =  (28M  )/l2 

'w  p 

is  either  a  solution  or  upper  bound. 

(2)  Same  as  mode  1  q  =  (28M  )/i2 

^  p 

Upper  bound  or  solution 


T)  Side  sway 

Upper  bound  or  solution 


q  =  (27.8  M  )/e 
P 


J)  Column  Mode  q  =  (44.4  M  )/  i2 

p 

Upper  boimd  or  solution 

(but  obviously  a  very  high  bound) 


\5)  Hinge  Mechanism 

(Degenerate  Mechanism)  No  energy  dissipation 


If  we  look  at  number  (S),  we  find 


q  =  (27.8  M  )/l2 
P 


and  this  is  either  an  upper  bound  or  a  solution.  Solving  for  the  equilibrium  required 

by  this  solution  we  find  an  overshoot  of  1.75  M  . 

P 


Hence  27.80/1.75  =  15.9  <  q  l2/ M  <  27.80 

P 

Lower  bound  Upper  bound 


This  very  first  approximation  boxes  the  answer  by  25%.  Going  through  several 
solutions,  we  find  the  solution  is  a  combination  of  Modes  2,  3,  and  5  plus  one-third  of 
Mode  4: 
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which  took  five  more  solutions  and  looked  promising.  But  adding  Mode  1 


produces 


q  =  (23.53  Mp)/£2 


and  a  10%  overshoot,  of  the  moment  diagram,  or 


23.53/1.10  =  21.40  <  q  Jl2/M  <  23.53 


and  finally,  success  is  achieved  with 


and 


which  is  little  different  from  the  median  of  the  first  attempt.  Hence,  the  bounds 
approach  a  very  strong  technique  of  estimating  the  correct  failure  load. 
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APPENDIX  B 


COMPUTER  ANALYSIS  FUNDAMENTALS 


COMPUTER  ANALYSIS  FUNDAMENTALS 


Computer  codes  designed  to  solve  static  analyses  of  structures  require  the 
idealization  of  the  two-  or  three-dimensional  structure  as  an  assemblage  of  finite 
elements  connected  at  joints.  The  technique  divides  a  structure  into  small  elements 
with  easily  defined  stress  and  deflection  characteristics.  Corners  and,  in  some 
instances,  midpoints  of  these  elements  are  the  locations  known  as  grid  points  or 
nodes  where  deflections  are  calculated  through  the  use  of  known  stiffness  properties, 
applied  loads,  and  boundary  conditions  or  displacement  constraints.  This  finite 
element  method  is  based  on  arrays  of  large  matrix  equations  that  appear 
complicated;  however,  the  method  is  based  on  fundamentally  simple  concepts 
involving  basic  stiffness  and  deflection  equations. 

A  model  must  be  constructed  that  accurately  represents  the  structural  stiffness 
and  displacement  constraints  and  provides  for  the  application  of  equivalent  loadings. 
The  input  data  are  processed  and  the  computer  output  provides  nodal  displacement, 
member  forces,  reactions,  and  stresses.  Most  large-scale  commercial  programs  also 
have  post-processing  routines  to  allow  the  selection  of  specific  output  data.  Stress 
contour  plots  can  be  generated  for  various  sections  through  the  model.  Deformed 
shape  plots  are  also  available  which  show  model  displacements  that  are  greatly 
exaggerated  for  clarity. 

The  computation  of  dynamic  structural  behavior  requires  the  solution  of  the 
dynamic  equations  of  equilibrium.  These  equations  may  be  written  in  the  following 
matrix  form: 


[M]{X}**+  [c]  {  X}**  +  [k]  (X)  *  {  F(t)  } 


(1) 


** 


X 


a2  x 

a  t2 


and  X 


9  X 

9  t 


where  t  represents  time. 
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where  [M]»  [C],  and  [K]  are  matrices  representing  the  mass,  damping,  and  stiffness 
of  the  structure  respectively.  The  vectors  {X},  {X},  and  {X}  represent  the 
accelerations,  velocities,  and  displacements  of  the  joint  degrees  of  freedom.  The 
vector  f  F(t)}  represents  the  dynamic,  time-varying  applied  forces.  Eq.  1  states 
that  the  sum  of  the  inertial,  damping,  and  elastic  resistance  forces  of  the  structure 
must  be  equal  to  the  applied  loads  at  each  point  in  time. 


Each  unrestrained  joint  degree  of  freedom  in  the  structure  is  represented  by  an 
equation  in  matrix  Eq.  1.  In  many  instances,  the  forces  acting  at  some  of  the 
degrees  of  freedom  are  small  and  do  not  significantly  influence  the  behavior  of  the 
structure.  The  equations  representing  these  degrees  of  freedom  may  be  eliminated 
in  order  to  reduce  the  total  degrees  of  freedom  of  the  structure  prior  to  solution. 
Two  techniques  are  used  for  reducing  the  dynamic  degrees  of  freedom.  These 
techniques  are  called  STATIC  and  GUY  AN  CONDENSATION.  To  understand  the 
difference  between  STATIC  and  GUYAN  CONDENSATION,  consider  the  undamped 
partitioned  form  of  the  Eq.  1. 
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The  subscripts  C  and  D  in  the  above  equation  denote  the  degrees  of  freedom  to  be 
condensed  and  the  dynamic  degrees  of  freedom  to  be  retained  for  solution 
respectively. 


Both  STATIC  and  GUYAN  CONDENSATION  assume  that  there  is  a  linear 
transformation  between  the  dynamic  degrees  of  freedom  and  the  condensed  degrees 
of  freedom: 

{\}  *  [T] 

{Xc}  -  [T] 


(3) 
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where  [T]  is  the  transformation  which  is  independent  of  time.  Using  Eq.  3,  both 
the  dynamic  and  condensed  displacements  and  accelerations  may  now  be  written  as  a 
function  of  only  the  dynamic  degrees  of  freedom: 


(4) 


and 


(5) 


where  [i]  is  the  identity  matrix. 

r  t'"' 

Substituting  Eq.  4  and  5  into  Eq.  2  and  premultiplying  by  [i]  [T]  ;to 
preserve  symmetry  yields: 
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Performing  the  indicated  operations  in  Eq.  6  results  in  the  following  two 
equations: 
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and 


M  ]  (X  }  +  m  P  {x  }  +  [k  1  .{x  }  +  “k  1  {x  }  =  {f  } 

CD’  ‘  DJ  cc:  1  C1  ;  CD]  1  DJ  CCl  1  CJ  1  C; 


(8) 


Using  Eq.  8  and  further  assuming  that  the  inertial  forces  do  not  affect  the 
relationship  between  {XQ}  and  { X^ }  (i.e.,  they  may  be  neglected),  Eq.  8  may  be 
solved  for  {  X_}  in  terms  of  {  X  },  neglecting  {  F  }  since  it  does  not  affect  the  free 
vibration  problem: 
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Thus, 


[T]  =  K 


(10) 


jyccj  cdj 

Substituting  Eq.  10  into  Eq.  7  results  in  the  following  condensed  equations  of  motion. 


[M]  {Xj  +  M  {X])}  -  [F] 
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(13) 
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Eq.  11  represents  the  condensed  equations  of  motion  for  both  STATIC  and  GUYAN 
CONDENSATION.  The  matrix  Ck]  in  Eq.  13  is  the  same  for  both  techniques. 
GTSTRUDL  restricts  the  applied  loadings  to  be  applied  only  to  the  dynamic  degrees 
of  freedom  (i.e..  {F_  }  =  0).  Therefore, 


{F}*  =  {fd}  (15) 

for  both  methods.  The  form  of  [  M]  is  different  in  the  two  techniques.  GUYAN 
CONDENSATION  uses  [M]  as  given  by  Eq.  12  while  STATIC  CONDENSATION 
utilizes  only  of  Eq.  12  for  [M]  .  A  physical  interpretation  of  the  difference 

in  the  two  techniques  is  that  GUYAN  CONDENSATION  redistributes  the  mass  of  the 
structure  including  the  mass  associated  with  the  condensed  degrees  of  freedom  while 
STATIC  CONDENSATION  considers  only  the  mass  associated  with  the  dynamic 
degrees  of  freedom.  GUYAN  CONDENSATION  is  the  recommended  technique 
unless  a  very  large  portion  of  the  mass  of  the  structure  is  associated  with  only  the 
dynamic  degrees  of  freedom.  The  only  advantage  of  STATIC  CONDENSATION  is 
reduced  computational  effort  required  to  form  [M]  . 


The  vibration  of  an  elastic  structure  is  assumed  to  be  harmonic  in  nature.  The 
total  vibration  of  an  elastic  structure  may  be  considered  to  be  a  linear  combination 
of  the  characteristic  harmonic  mode  shapes  denoted  by  {$} .  Thus,  the  displacement 
and  acceleration  of  any  characteristic  mode  of  vibration,  i,  may  be  expressed  in  the 
following  manner: 

{x}^  =  1$}^  sin  0)^  t  (16) 

{x}^  =  -{ <t>} ±  w*  sin  a)  t  (17) 
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where  oj.  is  the  frequency  of  vibration  in  the  i**1  mode.  The  undamped  equation  of 
free  vibration  is 

[k]  {x}  =  -[m]  {x}  (18) 

Substituting  Eqs.  16  and  17  into  18  yields, 

[Kl{$}i  =  uj  [M]  {$}±  (19) 

Eq.  19  is  therefore  an  eigenvalue  problem  in  nonstandard  form  in  which  the 
eigenvalues  are  the  squares  of  the  frequencies  and  the  eigenvectors  are  the 
vibrational  modes  associated  with  the  frequencies.  Considering  all  possible 
frequencies  and  mode  shapes,  Eq.  19  may  be  written  as 

[k]  [*]  -  [n]  [m]  O]  (20) 

where  [fi]  is  diagonal  containing  the  eigenvalues  on  the  diagonals  in  order  of 
increasing  magnitude.  The  matrix  [$]  contains  the  eigenvectors  stored  columnwise 
in  an  order  corresponding  to  the  eigenvalues  in  [fi].  The  matrix  [ft]  is  often  called 
the  spectral  matrix  while  [<&]  is  called  the  modal  matrix. 

Eq.  20  must  be  transformed  to  standard  eigenvalue  form  ([  D]  [r]  =  [fi]  [r])  in 
order  to  compute  the  eigenvalue  and  eigenvector  solution.  An  orthogonal  similarity 
transformation  is  used  for  the  transformation  in  order  to  preserve  symmetry  in  [D]. 
The  eigenvalues  of  the  original  equation  are  invariant  under  this  transformation 
while  the  eigenvectors  are  transformed.  For  lumped  mass,  the  transformation 
assumes  a  new  eigenvector  basis  such  that 

[*]  -  Cm]"1*  [r]  (21) 

where 

[M]  -  [M]*  [M]1*  (22) 
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The  matrix  [M]^  is  defined  to  be  the  square  root  of  all  elements  in[  M].  Since  [M] 
is  a  diagonal  matrix  for  lumped  mass,  [M]  is  the  square  root  of  only  the  diagonal 
elements.  The  matrix  C  M]  in  Eq.  21  is  a  diagonal  matrix  where  the  diagonal 
elements  are  one  over  the  square  root  of  the  diagonal  elements  in  [M]. 
Therefore, 


[r]  =  [Mlf5  [*] 


(23) 


Substituting  Eq.  21  into  20  yields 

[K]  [M]-55  [r]  =  [n]  [M]  [M]-1*  [r]  (24) 

Now,  substituting  Eq.  22  into  Eq.  24  produces 


[k]  Cr]  =  [n]  [mT5  [m]^  [m]'*5 

Multiplying  by  EmJ  and  realizing  that  [M]  1  [M] 
yields 


[r]  (25) 

*  is  equal  to  the  identity  matrix 


[d]  [r]  =  [n]  [r]  (26) 

where  [D]  is  symmetric  and  is  defined  as 

[D]  =  [M]'*5  [K]  W*  (27) 

Eq.  26  is  the  eigenvalue  problem  in  standard  form.  Again,  the  eigenvalues  in  [ft] 
have  remained  invariant  under  the  transformation  while  the  eigenvectors  have  been 
transformed  using  Eq.  23. 

If  the  response  of  a  structure  due  to  a  transient  loading  is  desired  as  a  function 
of  time,  a  time  history  analysis  must  be  performed.  For  the  time  history  analysis  a 
modal  superposition  method  is  used.  The  modal  superposition  method  requires  the 
solution  of  Eq.  19  for  the  frequencies  and  mode  shapes.  The  mode  shapes,  [$>],  can 
be  shown  to  be  orthogonal  with  respect  to  the  mass  and  stiffness  matrices.  The  mode 
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shapes  may  be  normalized  to  produce  the  following  results: 

[4]T  [m]  [4]  =  [l]  (28) 

When  normalized  according  to  Eq.  28,  the  mode  shapes  are  said  to  be  "normalized  to 
unit  mass".  When  the  mode  shapes  are  normalized  to  unit  mass,  the  following 
relation  is  also  valid: 


l>]T  [k]  [*]  =  [fl]  (29) 

Eqs.  28  and  29  permit  the  equations  of  motion  to  be  uncoupled  in  a  dynamic  analysis. 
Prior  to  uncoupling,  the  equation  of  motion  in  physical  (structure)  coordinates  is 
given  by 

[M]  {x}  +  [c]  {x}  +  [k]  {x}  =  {F(t)}  (30) 

The  accelerations,  velocities,  and  displacements  in  Eq.  30  are  transformed  to  a 
different  coordinate  system: 

{X}  =  [4]  {U} ;  {x}  =  [4]  {U};  and  {x}  =  [$]  {u}  (31) 

where  [4]  is  the  modal  matrix  normalized  to  unit  mass  as  presented  in  Eq.  28.  The 
vectors  { U },  {U  j,  and  {U )  in  Eq.  31  are  the  accelerations,  velocities  and 
displacements  respectively  of  the  structure  in  what  is  defined  as  the  normal 
coordinate  system. 

T 

Substituting  Eq.  31  into  Eq.  30  and  premultiplying  by  [$]  yields 
[*]T  Cm]  [4]  {u}  +  [4]T  [C]  [4]  {U}  +  [$]T  [K]  [4]  {u}  =  [4]T  {F(t)}  (32) 

Substituting  Eqs.  28  and  29  into  Eq.  32  yields 

{U}  +  [4]T  [c]  [4]  {0}  +  [n]  {U}  =  [4]T  {F(t)}  (33) 
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T 

Note  that  Eq.  33  will  be  completely  uncoupled  if  [$]  [c]  [$]  is  a  diagonal  matrix. 

The  inclusion  of  damping  in  Eq.  33  creates  some  difficulty.  There  are  few 
theoretical  means  for  determining  the  damping,  and  experimental  investigations  have 
produced  data  for  only  a  small  number  of  existing  structures.  The  analyst  can, 
however,  usually  assess  a  percentage  of  critical  damping  for  each  mode  of  interest. 
This  percentage  of  critical  damping  for  each  mode  i  is  called  the  damping  ratio 
for  mode  i.  This  simplification  allows  Eq.  33  to  be  completely  uncoupled,  and  a 
separate  equation  written  for  each  mode  i: 

U.  +  2n.w.D  +  c^U.  =  f.(t)  (34) 

l  ill  ii  l 

where 

th  T 

f.(t)  is  the  loading  at  the  i  normal  coordinate,  =  4> .  F(t), 

1  th  1 

n.  is  the  damping  ratio  for  the  i  mode,  and 

1  th 

oj.  is  the  frequency  associated  with  the  i  mode. 

Eq.  34  may  now  be  solved  since  it  is  a  linear  differential  equation.  GTSTRUDL  uses 
the  Duhamel  Integral  solution  which  is  exact  except  that  the  loading  is  assumed  to 
be  piecewise  linear  between  integration  points. 

The  modal  superposition  method  requires  that  the  eigenproblem  be  solved  first 
and,  if  damping  is  present,  damping  ratios  must  be  specified. 


ANALYSIS  PROCEDURE  FOR  A  BUILDING  STRUCTURE  UNDER  BLAST  LOADING 

The  object  of  this  computer-assisted  analysis  is  to  predict  when  failure  of  a 
multi-story  building  occurs  because  of  the  effects  of  a  time-varying  blast  load  and 
also  to  indicate  what  the  collapse  mechanism  of  the  structure  will  be.  Normally, 
the  first  step  in  any  usual  structural  analysis  would  be  the  application  of  a  static 
load  in  the  direction  of  the  blast  load.  However,  unless  a  correlation  can  be  found, 
based  on  previous  experience  or  other  criteria,  between  the  blast  load  and  the 
statically  applied  (pseudo-static)  load,  this  method  proves  to  be  highly 
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unsatisfactory.  In  order  to  best  predict  response  it  is  thus  necessary  to  apply  the 
blast  load  forcing  function  to  the  structure  as  a  truly  time -varying  load,  and 
determine  the  resulting  behavior  of  the  structure  as  a  time-varying  function. 


Dynamic  analyses  are  usually  conducted  within  the  range  of  elastic  behavior  of 
the  system  under  consideration.  A  collapsing  structure  is  outside  of  this  range  for 
the  obvious  reason  that  collapse  does  not  occur  in  the  elastic  domain.  Prior  to 
collapse  of  a  structure  several  behavioral  phases  must  be  satisfied.  Initially,  the 
structure  behaves  elastically  until  moments  and  forces  are  generated  in  structural 
components  to  cause  local  yielding,  i.e.,  localized  behavior  outside  the  elastic  range 
of  the  structural  material.  Plastic  hinges  are  formed  in  individual  members,  which 
will  now  resist  a  constant  moment  only,  independent  of  increased  bending  strain, 
except  for  a  small  increase  due  to  what  is  normally  referred  to  as  strain  hardening. 
If  a  sufficient  number  of  plastic  hinges  have  formed  to  cause  the  structure  to  behave 
as  a  mechanism,  collapse  is  imminent.  In  the  analysis  two  distinct  structural 
behavior  patterns  were  considered:  first,  elastic  behavior  until  a  sufficient  number 
of  columns  and  beams  had  reached  the  limits  of  their  elastic  range  to  cause  the 
structure  to  develop  a  failure  mechanism;  second,  plastic  behavior  until  total 
collapse  of  the  failure  mechanism  occurs  at  the  fracture  deformation  end  of  the 
plastic  range. 

The  plastic  behavior  was  modeled  using  a  reduced  modulus  of  elasticity  for 
those  structural  members  that  reached  the  limit  of  their  elastic  behavior.  External 
moments  equal  to  the  plastic  hinge  moments  were  applied  at  those  joints  where  a 
plastic  hinge  had  formed  as  part  of  the  failure  mechanism.  This  unique  method  of 
modeling  pseudo-plastic  behavior  retains  stability  of  the  structure,  even  though  a 
mechanism  has  formed,  and  also  provides  a  means  of  representing  the  strain 
hardening  of  the  structural  material.  To  retain  continuity  of  behavior  in  the 
transition  between  the  elastic  and  the  plastic  phase,  the  plastic  phase  was  started 
with  initial  conditions  of  displacement  and  velocity  equal  to  those  displacements  and 
velocities  at  the  end  of  the  elastic  phase,  see  Figure  B-l. 
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Fig.  B-l .  Computer  Model  for  Elasto-Plastic  Behavior. 


REPRESENTATION  OF  DUCTILITY  AND  ELEMENT  FAILURE  CRITERIA 


In  Section  4  -  Subtask  3b,  the  element  strength  and  recommended  ductility 
features  (ratio  of  failure  deformation  to  yield  deformation)  were  given  for  both  steel 
and  reinforced  concrete  structural  members.  In  order  to  incorporate  these 
properties  into  the  computer  model,  it  is  necessary  to  express  the  failure  criteria  in 
terms  of  bending  moments  rather  than  member  section  deformations  (rotational 
angles).  This  is  because  the  basic  computer  output  is  in  terms  of  element  forces 
(moments)  rather  than  deformations. 

COMPUTED  MOMENT  M 


H2TATICKI  DeF%^IWST)d>N  G 


For  example,  if  member  failure  is  defined  by  a  moment  of  1.5  M  in  the 

p 

inelastic  (1/10  El)  portion  of  member  deformation,  then  this  corresponds  to  a 
ductility  given  by 

0.5  Mp  =  0.5  El  fy  =  1/10  El  (u  0  ) 
or 

v  =  5 
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APPENDIX  C 

DEBRIS  GENERATED  LOADS 


DEBRIS  GENERATED  LOADS 


The  ability  to  predict  debris  generated  loads  depends  greatly  upon  the 
predicted  strength  of  exterior  and  interior  walls*.  According  to  this  reference, 
debris  prediction  and  the  resulting  loads  are  only  valid  for  accurate  predictions  of 
wall  and  partition  strength.  Data  concerning  building  response  to  blast  may  be 
qualitatively  useful  when  classified  in  terms  of  light,  moderate,  and  severe  damage, 
but  do  not  reflect  any  information  on  debris  production  or  distribution.  Thus,  data 
from  shock  tunnel  tests  and  the  Nevada  weapons  tests  were  used  to  construct  debris 
charts.  These  charts  predict  the  percent  of  building  material  that  will  become 
debris  at  various  incident  overpressures  for  a  specific  type  of  construction. 

The  intended  use  of  these  charts  is  to  provide  additional  data  to  develop  floor 
loadings  that  will  affect  the  frame  response.  In  the  case  of  complete  collapse,  the 
frame  response,  combined  with  debris,  will  load  the  roof  of  the  basement  shelter. 
For  any  event  less  than  that  to  cause  total  collapse,  the  following  general  scenario 
applies  for  a  multi-story  steel  or  reinforced-concrete  frame  structure,  with  many 
frangible  walls,  having  the  same  general  proportions  as  the  Peachtree  Building. 

The  exterior  wall  facing  the  blast  will,  upon  failure,  transfer  its  horizontal  and 
vertical  energy  to  the  floor  below  through  friction  and/or  collision  with  interior 
partitions  and  objects.  The  exterior  walls  on  the  back  side,  if  failed,  will  fall  onto 
an  adjacent  site.  At  40  psi  incident  overpressure,  and  for  a  blast  of  long  duration, 
considerable  debris  from  interior  and  exterior  walls  may  leave  the  site  as  well. 
However,  up  to  a  point,  long  duration  blast  waves  will  also  greatly  increase  the  drag 
loading  —  thereby  increasing  the  structural  damage  and  the  production  of  debris. 
Therefore,  considering  total  wall  weight  as  debris  for  more  than  3  psi  overpressure  is 
conservative. 


*  Edmunds,  James  E.,  Structural  Debris  and  Building  Damage  Prediction  Methods, 

URS  686-5,  URS  Research  Company,  Burlingame,  CA,  June  1968. 
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The  vertical  energy  imparted  to  the  floor  below  can  be  expressed  as  the  weight 
of  the  element  times  the  height  through  which  it  falls.  However,  the  resulting  force 
is  very  small  compared  to  existing  vertical  forces,  especially  when  the  weight  of  the 
exterior  wall  is  considered  to  be  approximately  evenly  distributed  over  the  floor 
area.  The  horizontal  force  vector  generated  by  the  debris  can  be  determined  from 
the  drag  loading  expressed  as 


F(t)  =  Cd  q(t)  A 

The  dynamic  pressure  pulse,  q(t),  acts  on  the  exposed  debris  in  the  same 
manner  as  on  the  exposed  frame  with  a  area  A  and  a  drag  coefficient  C^.  If  the 
coefficient  of  friction,  y,  times  the  debris  weight  exceeds  F(t)  above,  the  debris  will 
not  move.  As  long  as  F(t)  is  greater  than  the  debris  weight  times  y,  the  debris  will 
slide  across  the  floor  generating  a  horizontal  force  on  the  floor,  generally  until  it  is 
carried  out.  Assuming  y  is  independent  of  velocity  (i.e.,  nearly  constant),  the 
maximum  horizontal  force  generated  by  debris  will  be 

F  =  weight  *y 
max  &  s 

for  a  duration  less  than  t+.  Table  C-l  shows  the  approximate  values  of  yg  for 
different  surfaces. 


TABLE  C-l:  STATIC  COEFFICIENTS  OF  FRICTION  OF  DRY  MATERIALS 


Material 

»*s 

Concrete  to  concrete 

0.8 

Concrete  to  steel 

0.4 

Debris  on  concrete 

0.6 

The  procedure  to  determine  the  debris  generated  frame  loading  is  to  first  find 
the  incident  overpressure  and  duration  for  a  given  weapon  (see  Section  4,  Subtask 
3A)  then,  determine  weight  of  wall  material  assuming  that  it  becomes  debris.  Use 
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Table  C-2  for  a  given  building  type  and  determine  the  volume  of  structural  materials. 

This  volume  times  the  unit  weight  gives  the  weight  of  the  structural  debris.  The 

contents  (interior  partitions,  etc.)  can  also  contribute  significantly.  Therefore, 

determine  volume  of  contents  from  the  type  of  occupancy  in  Table  C-3  and  multiply 
3 

by  60  lbs  per  ft  to  obtain  contents  weight,  and  add  this  to  the  weight  of  structural 
debris  for  total  debris  weight.  Multiply  total  weight  (walls  and  contents)  by  the 
appropriate  coefficient(s)  of  friction  (from  Table  C-l)  and  add  this  to  the  structural 
frame.  Mathematically  this  procedure  is  expressed  as: 


Frame  Loading  = 


V  •  w(wt/ft^)  +  KA  u 
c  p  s 


where 


Vc  is  volume  of  structural  content 
w  is  the  weight  per  unit  volume  of  structural  material 
K  is  a  contents  coefficient  from  Table  C-3 
Ap  is  the  plan  area,  and 
ys  is  coefficient  of  friction  from  Table  C-l. 


Of  course,  this  expression  is  for  one  story;  the  same  procedure  can  be  used  for 
each  story  in  a  multi-story  building. 


DEBRIS  LOADING  DATA 

In  Subtask  3A  (see  Section  4)  available  research  establishes  the  failure  of 
exterior  walls  due  to  blast  pressures  at  the  range  of  5  to  15  psi.  In  the  computer 
analysis  this  100%  debris  weight  is  represented  by  an  increase  in  the  assigned  gravity 
load  (or  dead  load)  on  the  structure  model. 

The  corresponding  friction  drag  force  of  the  moving  debris  action  on  the 

building  frame  and  floor  slabs  is  represented  by  the  use  of  a  high  "conservative"  drag 

coefficient  C .  in  the  evaluation  of  the  drag  pressure  portion  of  the  blast  loading 
d 

function. 


C-3 


TABLE  C-2 


STRUCTUSE  VOLUME  VS  BUILDING  TYPE 


BUILDING  TYPE 

VOLUME  FORMULA* 

1.  Wood  Frame  Residential 

V  =  A  •  height 

a.  1st  floor  slab  on 

ground 

[0.55  +  (N-l)  (0.525)1  A 

b.  1st  floor  on  std. 

P 

joists 

[0.7  +  (N-l)  (0.525)]  Ap 

2.  Steel  Frame  Industrial 

a.  Light  W/CI  sheathing 

0.02  A 

W/CA  sheathing 

0.087  AP 

b.  Heavy  W/CI  sheathing 

0.037  AP 

W/CA  sheathing 

0.095  AP 

P 

3.  Load-Bearing  Masonry  with 

or  without  reinforcing  - 

0.12  V 

Combustible  interior 

framing 

4.  Reinforced  Concrete 

Shear- Wall 

a.  W/lt.  interior  panels 

0.07  V 

b.  W/ masonry  interior 

C 

panels 

0.12  V 

c 

5.  Multistory  Steel  and  Re- 

inforced  Concrete  Frame 

with  Earthquake  Design 

a.  W/lt.  interior  panels 

0.07  V 

b.  W/ masonry  interior 

panels 

0.11  V 

c 

6.  Multistory  Steel  and  Re- 

inforced  Concrete  Frame 

(non-earthquake  design) 

a.  W/lt.  interior  panels 

0.063  V 

b.  W/ masonry  interior 

c 

panels 

0.10  V 
c 

These  formulae  reflect  solid  volume  of  material  (i.e.,  no  void-void  ratio 
=  0) .  The  void  ratio  (usually  taken  as  unity)  is  best  applied  after  summa¬ 
tion  of  contributory  volumes.  This  minimizes  the  number  of  calculations  re 
quired  for  making  debris  volume  or  debris  depth  estimates. 


TABLE  C-2  (contd) 


STRUCTURE  VOLUME  VS  BUILDING  TYPE 


BUILDING  TYPE 

VOLUME  FORMULA 

7.  Light  Reinforced  Concrete 

Shear-Wall  (single  story) 

a.  Concrete  roof  w/lt. 

interior  panels 

0.07  V 

b.  Concrete  roof  w/ma- 

C 

sonry  interior  pan- 

els 

0.075  V 

c.  Mill  roof  w/lt.  int. 

C 

panels 

0.037  V 

d.  Mill  roof  w/masonry 

interior  panels 

0.05  V 

c 

LEGEND: 

V  contained  volume 
c 

A  plan  area 
P 

S&P  stucco  exterior  plaster  interior 
W&P  wood  exterior  plaster  interior 
W  all  wood 
Cl  corrugated  iron 
CA  corrugated  asbestos 


TABLE  C-3 


BUILDING  CONTENTS  LOADS  AND  VOLUME  FACTORS 


OCCUPANCY 

PSF 

TOTAL 

VOLUME  FACTOR  K 
(v '  “p  )* 

TOTAL 

1  . 

Apts,  and  Residential 

5 

0.625 

2. 

Auditoriums  and  Churches 

1.5 

0.25 

3. 

Ga  ra  ge 

a  .  S  to  ra  ge 

15 

0.75 

b.  Repair 

11 

0.55 

4. 

Gymnasium 

0.5 

0.09 

5. 

Hospitals 

3 

0.375 

6  . 

Hotels 

5 

0.625 

7. 

Libraries 

26 

0.75 

8. 

Manufacturing 

a.  Comb.  Mdse,  fabrics, 

furniture 

18 

1.8 

b.  Incombustible 

11 

0.55 

9. 

Offices 

12 

1  .2 

10. 

Printing  Plant 

a.  Newspaper 

23 

0.9 

b .  Books 

60 

1.7 

11  . 

Schools 

11 

1.6 

12. 

Storage 

a .  Gen .  Mdse . 

35 

6 

b.  Special 

** 

13. 

Stores 

a.  Retail  Dept. 

12 

2 

b.  Wholesale 

16 

2.7 

14. 

Restaurant 

3.5 

0.6 

V  =  Volume  in  cubic  feet  *  KA 

P 

=  Plan  area  in  square  feet 
P 


06 


It  is  most  appropriate  here  to  note  a  lack  of  research  information.  Specif¬ 
ically,  all  of  the  exterior  wall  blast  research  is  in  the  range  of  1  psi  <  pQ<  10  psi 
loading.  The  area  of  interest  in  this  research  is  in  the  range  of  30  psi_£  p^£  50  psi, 
and  it  is  most  difficult  to  extrapolate  and  forecast  wall-frame  interaction  forces  and 
the  resulting  debris  drag  load  effects  in  this  high  pressure  range.  There  is  certainly 
a  need  for  both  analysis  and  testing  research  to  bridge  this  knowledge  gap  for  high 
blast  pressure  effects  on  debris  and  debris-structure  interaction. 
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